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Abstract 

We study Brownian motion and the associated Langevin equation in AdS/CFT. The 
Brownian particle is reahzed in the bulk spacetime as a probe fundamental string in an 
asymptotically AdS black hole background, stretching between the AdS boundary and 
the horizon. The modes on the string are excited by the thermal black hole environment 
and consequently the string endpoint at the boundary undergoes an erratic motion, 
which is identified with an external quark in the boundary CFT exhibiting Brownian 
motion. Semiclassically, the modes on the string are thermally excited due to Hawking 
radiation, which translates into the random force appearing in the boundary Langevin 
equation, while the friction in the Langevin equation corresponds to the excitation on 
the string being absorbed by the black hole. We give a bulk proof of the fluctuation- 
dissipation theorem relating the random force and friction. This work can be regarded 
as a step toward understanding the quantum microphysics underlying the fluid-gravity 
correspondence. We also initiate a study of the properties of the effective membrane or 
stretched horizon picture of black holes using our bulk description of Brownian motion. 
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1 Introduction 



One of the interesting problems in statistical mechanics concerns the understanding of the 
origin of macroscopic dissipation and the approach to thermal equilibrium from microscopical 
point of view. Conventionally, given a statistical system in the thermodynamic or hydrody- 
namic limit, we imagine the collisions between the microscopic constituents of our system as 
being responsible for both of these macroscopic phenomena. This kinetic theory based pic- 
ture is firmly anchored on the basic idea of Brownian motion — in 1827, the botanist Robert 
Brown observed |T]] under a microscope that tiny pollen particles suspended in water undergo 
incessant irregular motion, which became known as the Brownian motion.]^ As is well-known 
now, this peculiar motion is due to collisions with the fluid particles in random thermal mo- 
tion. Therefore, any particle immersed in fluid at finite temperature exhibits such Brownian 
motion, from a small pendulum suspended in a dilute gas to a heavy particle in quark- 
gluon plasma. This universal phenomenon suggests that the interaction with microscopic 
constituents is responsible for dissipation and thermalization on macroscopic scales. 

Since its advent, the holographic AdS/CFT correspondence 0,|,||,[T^ has been exploited 
to study the physics of non-Abelian quark-gluon plasmas at finite temperature from bulk 
gravitational physics, and vice versa. The dual gravitational description of strongly coupled 
gauge theories provides an efficient way to study the thermodynamic properties and the phase 
structure of the gauge theory. More recently, it has become clear that one can also exploit the 
gravitational description to understand the hydrodynamic regime of the quark-gluon plasma. 



as was originally proposed in |TT[| and has been significantly developed afterwards (see [|T2[ 
and references therein for earlier work on hydrodynamics in the AdS/CFT context). Namely, 
the long-wavelength physics described by a hydrodynamical Navier-Stokes equation on the 
boundary side is holographically dual to the long-wavelength fiuctuation of the horizons of 
asymptotically AdS black hole spacetimes on the gravitational side. This correspondence 
allows for a detailed quantitative study of the plasma from the bulk, and vice versa. It is thus a 
natural question to ask whether one can obtain a holographic description of Brownian motion, 
which is one step towards the microphysics underlying thermodynamics and hydrodynamics. 
The aim of this paper is to answer this question in the affirmative. 

One intrinsic reason to be interested in Brownian motion within a holographic setting is 
to better understand the microscopic origin of the thermodynamic properties of black holes. 
It has become clear from the formulation of the AdS/CFT correspondence that one has an 
in-principle solution to the problem of quantum dynamics of black holes: we only need to solve 
the problem phrased in terms of the dual field theory variables. However, it is fair to say that a 
concrete quantitative understanding of the physics in these contexts is still lacking. One of the 
most useful playgrounds for understanding the quantum behavior of black holes has been the 



^Classic reviews on Brownian motion are [^,^,^. For a more complete list of references, see e.g. 
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arena of supersymmetric black holes |13, 3- Here we not only understand in many cases the 
microscopic origin of black hole entropy, but also in a number of cases have a bulk picture of 
the nature of the microscopic states making up the black hole degeneracy. In fact, from these 
various analyses, there emerges a rather intriguing picture of a quantum black hole — the black 
hole microstates form a sort of spacetime foam that replaces the region inside the horizon. 
Any single microstate is horizon-free, but the typical microstates are expected to exhibit 
the characteristic features of black hole spacetimes, which has been confirmed explicitly for 
some concrete systems; see [15,|16|, 1^ 1^] for reviews. Given this state of affairs one might 
probe these microstates beyond equilibrium thermodynamics and ask how the ensemble of 
them leads to dissipation and thermalization seen in a thermal medium. Understanding the 
description of Brownian motion seems then a natural step towards getting a handle on the 
problem. 

Conversely, as mentioned above, the AdS/CFT correspondence has been immensely use- 
ful in understanding many qualitative (and sometimes quantitative) features of quark-gluon 
plasmas. The famous lower bound on the ratio of shear-viscosity to entropy density for rela- 
tivistic hydrodynamic systems, rj/s > l/Air (see also a review [^), has certainly played 
an important role in obtaining a quantitative understanding of the dynamics of the quark- 
gluon plasma produced at RHIC. Furthermore, studies of the motion of quarks, mesons, 
and baryons in the quark-gluon plasma have been carried out in the holographic framework 

,^^,|2^, by considering the dynamics 



starting with the seminal papers pO| , pi| 



of probe strings and D-branes in asymptotically AdS black hole spacetime — for a sample of 



recent reviews on the subject, see [|28[. The general philosophy in these discussions was to use 
the probe dynamics to extract the rates of energy loss and transverse momentum broadening 
in the medium, which bear direct relevance to the physical problem of motion of quarks and 
mesons in the quark-gluon plasma. In such computations, the motion of an external quark in 
the quark-gluon plasma is assumed to be described by a relativistic Langevin equation [ p9| . 
In the most basic form, the Langevin equation is parametrized by two constants: the fric- 
tion (drag force) coefficient 7 and the magnitude of the random force k.0 Furthermore, the 
random force is assumed to be white noise. By using the AdS/CFT realization of external 
quarks, Refs. |2^,|2^ determined the friction coefficient 7, while Refs. |l2^,|25|,|27|l computed 
the random force k.| 

Therefore, one can say that the most basic data of the Langevin equation describing 



^ More precisely, in the relativistic case, the random force has different magnitudes kl and kt in the 
directions transverse and longitudinal to the momentum p. In the non-relativistic limit p ^ 0, they are equal: 
hl = Kt- The parameters 7 and kl are related to each other by the Einstein relation, under the assumption 
that the Langevin dynamics holds and gives the Jiittner distribution e~P^ . On the other hand, k,t is an 
independent parameter jH, 25 1 . 

■^More precisely, [^ , |22|] computed 7 and |25, 2^ computed kt, both in the relativistic case (the computation 
of was nonrclativistic) . The longitudinal component does not have to be computed independently, 
since it is related to 7 by the Einstein relation. See also footnote |[ 
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Brownian motion in the CFT plasma are already available. However, rather than taking such 
approaches which are phenomenological in some sense, one could study more fundamental 
aspects of Brownian motion in the AdS/CFT context. For example, in the first place, why 
does an external quark exhibit Brownian motion, and why is the motion described by a 
Langevin equation? While the domain of validity of the Langevin equation is clear from 
the previous results on the drag force, can we identify the origins of the Brownian motion 
approximation from the bulk gravitational description? What is the bulk meaning of the 
fluctuation-dissipation theorem relating 7 and k? The main purpose of the current paper is 
to elucidate the AdS/CFT physics of Brownian motion, by addressing such questions. For 
example, the computation of the random force in [^,^,^ using the GKPW prescription 
does not explain what the bulk counterpart of the random force is. We will see that it 
corresponds to a version of Hawking radiation in the bulk.^ 

Since we want to model Brownian motion, we need a gravitational analog of a particle 
immersed in a thermal medium. In the boundary field theory a natural particle is a test 
quark of large but finite mass immersed in the quark-gluon plasma. This is realized in 
the dual gravitational picture by introducing a fundamental string in the Schwarzschild-AdS 
background. The endpoint of the string at the boundary then corresponds to the test quark 
which undergoes Brownian motion; see Figure [l|. 



boundary 






string 



horizon ^ 

i ^ .r=rr 



Figure 1: The bulk dual of a Brownian particle: a fundamental string hanging 
from the boundary of the AdS space and dipping into the horizon. The AdS 
black hole environment excites the modes on the string and, as a result, the string 
endpoint at infinity moves randomly, corresponding to the Brownian motion on 
the boundary. 

We will use this simple picture of a probe fundamental string in a black hole background 
to "derive" the Brownian motion which the string endpoint on the boundary undergoes.^ 

^For an earlier discussion of the relation between Hawking radiation and diffusion the context of AdS/QCD, 
see 1^. 

preliminary discussion of the fluctuations of a fundamental string in an asymptotically AdS black hole 
background can be found in ]3l|] . 



4 



The basic idea is to quantize the fluctuations of the string world-sheet about a classical 
solution, which in the situations of interest corresponds to a straight string hanging down 
from the boundary. Since the bulk geometry has an event horizon, the induced metric on 
the string world-sheet also corresponds to a black hole geometry and the problem of studying 
fluctuations reduces to the dynamics of two dimensional quantum fields in curved spacetime. 
By quantizing the fluctuations we relate the quantum modes of the string to the boundary 
endpoint. This mapping in principle allows one to use the correlation functions of the position 
of the string endpoint on the boundary to recover the excitation spectrum of the string 
world-sheet. Assuming the validity of the semi classical approximation, we can then relate the 
thermal physics of the Hawking radiation to the Brownian motion of the string endpoint and 
derive the Langevin equation for the boundary dynamics. 

The organization of the rest of the paper is as follows. In section we set the stage for 
our discussion by first reviewing the Langevin equation describing Brownian motion in the 
field theory context. We then turn to a holographic realization of Brownian motion in terms 
of the dynamics of a probe fundamental string stretching from the boundary to the horizon 
of an asymptotically AdS black hole geometry in d dimensions. We write down the explicit 
relation between the boundary and bulk quantities associated with the holographic Brownian 
motion. This boundary-bulk relation can be explicitly worked out at the semiclassical level, 
which we turn to in section |^, focussing on the simple case of three dimensional spacetimes. 
There, we assume that modes on the string are thermally excited due to Hawking radiation 
and we derive the Langevin dynamics exhibited by the boundary Brownian particle. The 
friction and the random force appearing in the Langevin equation are related to each other 
by the fluctuation-dissipation theorem. In section ^, we study this theorem from the bulk 
viewpoint and give a bulk proof of it in the general case. In section |^, we generalize the 
discussion in section ^ for = 3 to general dimensions. Despite being unable to quantize the 
modes on the string analytically in this case we nevertheless show that at small frequencies 
we recover the Langevin equation. In section ^, we study whether the bulk Brownian motion 
of the fundamental string can be interpreted as being caused by a suitable movement of the 
string endpoint on the horizon, the idea being that the endpoint is randomly excited by 
the stringy gas living on a membrane just outside the horizon, much as in the spirit of the 
membrane paradigm. We also provide a preliminary discussion of how to use our setup to 
study microscopic properties of the stretched horizon. Ultimately, we would like to directly 
probe properties of the quasi-particles that make up the stringy gas living at the stretched 
horizon, but that is beyond the scope of the present paper. Section]^ is devoted to a discussion. 
Some of the relevant technical details are collected in the Appendices. 
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2 Holographic Brownian motion 



To set the stage for our discussion we begin with a brief review of the Langevin dynamics that 
describes the Brownian motion. This discussion will be the field theoretic, or boundary, side 
of the story in the AdS/CFT context. Turning to the corresponding bulk description, we will 
then describe how one can set up the problem of studying the motion of a Brownian particle 
in a thermal medium in terms of a probe string in an asymptotically Schwarzschild-AdS black 
hole background. 

2.1 Brownian motion and Langevin dynamics 

Let us begin with the Langevin equation, which is the simplest model describing a non- 
relativistic Brownian particle of mass m in one spatial dimension: 

pit) = --fop{t) + R{t), (2.1) 

where p = mx is the (non-relativistic) momentum of the Brownian particle at position x, and 
' = d/dt. The two terms on the right hand side of ( |2.1[ ) correspond to friction and a random 
force, respectively, and 70 is a constant called the friction coefficient. One can think of the 
particle as losing energy to the medium due to the friction term and concurrently getting a 
random kick from the thermal bath modeled by the random force, which we assume to be 
white noise with the following average: 

{R{t)) = 0, {R{t)R{t')) = 6{t - t'), (2.2) 

where kq is a constant. The separation of the force into frictional and random parts on the 
right hand side of (|2.1|) is merely a phenomenological simplification — microscopically, the two 
forces have the same origin (collision with the fiuid constituents). 

Assuming equipartition of energy, (mx^) = T, with T the temperature,^ one can derive 
the following time evolution for the displacement squared p[: 



(sit)') ^ {[x{t) - x(0)]2) = ^ (7ot - 1 + e-^°*) ^ <J 

7o 



T 2 




1 


— t^ 


(t < 




m 




7o 


2Dt 


(t > 


1 






7o 



(2.3) 



where the diffusion constant D is related to the friction coefficient 70 by the Sutherland- 
Einstein relation, 

D = (2.4) 

7om 



^We shall work in units where the Boltzmann constant fc^ = 1. 
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We can see that in the balhstic regime, t <^ I/70, the particle moves inertially (s ~ t) with 
the velocity determined by equipartition, x ~ ^/T/m, while in the diffusive regime, t ^ I/70, 
the particle undergoes a random walk (s ~ ^/t). This is because the Brownian particle must 
be hit by a certain number of fluid particles to get substantially diverted from the direction 
of its initial velocity. The crossover time between the two regimes is the relaxation time 

1 , , 

^relax ~ , (2.5) 

7o 

which characterizes the time scale for the Brownian particle to forget its initial velocity and 
thermalize. One can also derive the relation between the friction coefficient 70 and the size of 
the random force Kq 

70 = (2.6) 

which is the simplest example of the fluctuation-dissipation theorem and arises due to the 
fact that the frictional and random forces are of the same origin. 

In n spatial dimensions, p and R in ( p.l|) are generalized to n component vectors and (|2.2| ) 
is generalized to 

{R,{t)) = 0, {R,{t)R,{t')) = 5,, 5{t - t'), (2.7) 

where i, j = 1, . . . ,n. In the diffusive regime, the displacement squared goes as (s(t)^) ^ 
2nDt. The Sutherland-Einstein relation (|2.4| ) and the fluctuation-dissipation relation (|2.6| ) 
are independent of n. 

Now let us go back to the case with one spatial dimension (n = 1). The Langevin equation 
( p.l|) , ( |2.2|) captures certain essential features of physics, but nevertheless is too simple, for 
two reasons. It assumes that the friction is instantaneous and that there is no correlation 
between random forces at different times (Eq. ( |2.2| )). If the Brownian particle is not infinitely 
more massive than the fluid particles, these assumptions are no longer valid; friction will 
depend on the past history of the particle, and random forces at different times will not be 
fully independent. We can incorporate these effects by generalizing the simplest Langevin 



equation ( pTLf ) to the so-called generalized Langevin equation p2|,P3 



p{t) = - I dt'-iit - t')p{t') + R{t) + K{t). (2.8) 
J —00 

Now the friction term depends on the past trajectory via the memory kernel 'j{t). The random 
force is taken to satisfy 

{R{t)) = 0, {R{t)R{t')) = Kit - 1'), (2.9) 

where K{t) is some function. We have also now introduced an external force K{t) that can be 
applied to the system. 
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To analyze the physical content of the generalized Langevin equation we Fourier transform 
( p.8|) to obtain 

where p{uj), R{uj), K{uj) are Fourier transforms, e.g., 

/oo 
dtp{t)e'^\ (2.11) 
-oo 

while 7[cij] is the Fourier-Laplace transform: 

POO 

7[cj] = / dt-f{t)e''^K (2.12) 







If we take the statistical average of ( p.lO|) , the random force vanishes because of the first 
equation in ( |2.9| ), and we obtain 

{piu)) = ii{u)K{Lo), ^ , I . . (2.13) 

^{u}) is called the admittance. So, we can determine the admittance /i(co'), and thereby ^[uj], 
by measuring the response {p{uj)) to an external force. In particular, if the external force is 

K{t) = Koe-''^\ (2.14) 

then {pit)) is simply 

(p(t))=M^^)^oe-*"*. (2.15) 
For a quantity 0(t), define the power spectrum Io{^) by 

POD 

Ioiuj)= / dt{Oito)0{to + t))e'^\ (2.16) 



Note that {0{to)0(to + 1)) is independent of to in a stationary system. The knowledge of 
power spectrum is the same as that of 2-point function, because of the Wiener-Khint chine 
theorem 

{0{u)0{u')) = 27i6{uj + uj')Io{uj). (2.17) 
Now consider the case without an external force, i.e., K = 0. In this case, from (|2.10|) , 



pM = ^^M_, (2,18) 



Therefore, the power spectrum of p and that for R are related as 



/p(-) = TTT^r (2-19) 

|7 — %UJ\^ 
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Combining ( p.l5|) and ( p.l9|) , one can determine both 7(t) and K(t) appearing in the 
Langevin equation ( p.8|) and ( p.9|) separately. However, as we will discuss in section ^ these 
two quantities are not independent but are related to each other by the fluctuation-dissipation 
theorem, which is the generalization of the relation (|2.6|). 



For the generalized Langevin equation, what corresponds to the relaxation time (2.5) is 

1 



t 



relax 



dt-f{t) 



-f[uj = 0] 



fi{uj = 0). (2.20) 



If 7(t) is sharply peaked around t = 0, we can ignore the retarded effect of the friction term 
in (\2-8\) and write 

dt'^{t-t')p{t')^ I dt' ^{t') ■ Pit) = Pit). (2.21) 

J-rclax 



JO 



Then the Langevin equation reduces to the simple Langevin equation (|2.1| ) and it is clear that 
^reiax correspouds to the thermalization time for the Brownian particle. 

Another physically relevant time scale, the microscopic (or collision duration) time tcou, is 
defined to be the width of the random force correlator function K-it). Specifically, let us define 



CcoU 







rft^. (2.22) 

/€(0) 



If nit) = K(0)e~*/*™", the right hand side of this precisely gives tcou- This tcou characterizes 
the time scale over which the random force is correlated, and thus can be thought of as the 
time elapsed in a single process of scattering. In many cases, 

^relax ^ ^coU- (2.23) 

Typical examples for which ( p.23| ) holds are the case where the particle is scattered occasion- 
ally by dilute scatterers, and the case where a heavy particle is hit frequently by much smaller 
particles |^ . As we will see later, for the Brownian motion dual to AdS black holes, the field 
theories are strongly coupled CFTs and (|2.23|) does not necessarily hold. 



There is also a third natural time scale tmfp given by the typical time elapsed between two 
collisions. In the kinetic theory, this mean free path time is typically tcou ^mfp ^ ^reiax! 
however in the case of present interest, this separation no longer holds, as we will see. 

2.2 Bulk counterpart of Brownian motion 

The AdS/CFT correspondence states that string theory in AdS^ is dual to a CFT in (d — 1) 
dimensions. In particular, the (planar) Schwarzschild-AdS black hole with metric 

ds^, = C \-hir) de + dXl^ + dr\ /i(r) = 1 - ^JLY'^ (2.24) 



£2 
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is dual to a CFT at a temperature equal to the Hawking temperature of the black hole 

/? 47r£2 



1 (d-l)rH , , 

T = -= ^ ^ ;„ ^ . (2.25) 



In the above, i is the AdS radius, and t, Xct_2 = {X^, . . . ,X'^^^) G Mf^'^ are the boundary 
coordinates. 

In this black hole geometry (|2.24 ), let us consider a fundamental string suspended from 



the boundary at r = oo, straight down along the r direction, into the horizon ai r = th] see 
Figure |T]. In the boundary CFT, this corresponds to having a very heavy external charged 
particle. The Xa_2 coordinates of the string at r = oo in the bulk give the boundary position 
of the external particle. As we discussed above, such an external particle at finite temperature 
T is expected to undergo Brownian motion. The dual statement must be that the black hole 
environment in the bulk excites the modes on the string and, as the result, the endpoint of the 
string at r = oo exhibits a Brownian motion which can be modeled by a Langevin equation. 

We study this motion of a string in the probe approximation where we ignore its backre- 
action on the background geometry. We assume that there is no i?-field in the background, 
which is the case for AdSa based on Dl- and D5-branes, and AdSs based on D3-branes. If we 
take the string coupling Qs to be very small, the interaction of the string with the thermal gas 
of closed strings in the bulk of the AdS space can be ignored; the only possible region with 
appreciable interaction is near the black hole horizon which the string is dipping into. 



Let us slightly generalize ( |2.24| ) for a little while and consider the following metric: 



ds"^ = gf,uix) dx^ dx" + Gij{x) dX^ dX\ (2.26) 

where x'^ = t,r and I, J = 1, . . . ,d — 2. For the spacetimes of interest, both g^i, and Gjj 
are independent of X^ ^ Now, we stretch a string along the r direction and consider small 
fluctuation of it in the transverse directions X^ . The action for the string is simply the Nambu- 
Goto action in the absence of 5-field. In the gauge where the world-sheet coordinates are 
identified with the spacetime coordinates x^ = t, r, the transverse fluctuations X^ become 
functions of x^: X^ = X\x). If we expand the Nambu-Goto action up to quadratic order in 
X^ , we obtain 



5'ng = -77—7 / d'^x^/^det^ 
ZTra' J 



^ ^ d'xy^)g^'^{x)Guix)^^^S^^l (2.27) 



Ana' J dxt^ dx" 

where 7^,^ is the induced metric, g^'^ is the inverse of g^v, and g = det gfj^^- In the last line we 
dropped the constant term that does not depend on X^ . The quadratic approximation is of 



''We will mainly focus on planar black holes in AdS corresponding to thermal field theories on M"^ ^'^ when 
the transverse directions to the string are indeed Killing directions in the bulk. 
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course valid as long as the scalars do not fluctuate too far from their equilibrium value 
(taken here to be = O).0 In fact, this quadratic fluctuation Lagrangian for the world- 
sheet scalars ( |2.27 ) can be thought of as taking the non-relativistic limit; the Nambu-Goto 
action is after all a non-polynomial action in the velocities dtX^ and we are expanding in the 
regime |i9tX^| ^ 1. Therefore, we expect (and will see) that the dual Langevin dynamics on 
the boundary will also be a non-relativistic one, which is precisely what we reviewed in the 

(2) 

previous subsection. For most of the paper, we will use this quadratic action S^q to study 
the fluctuations of the string.]^ The equation of motion derived from ( |2.27| ) is 

1 



= V^iGjjd^X') 



-9 



d^i^g^'Gud^X- 



2.28 



where is the covariant derivative with respect to g^^. Note that this is not the same as 
the Klein-Gordon equation in the spacetime ( p. 261 ), which would involve not just 9^ = d/dx^ 
but also di = d/dX^ . 

Returning to the AdS black hole metric (|2.24|) , we focus first on the AdSs {d = 3) case for 
simplicity (we will discuss AdS^ with general d in section |^) and study the motion of a string 
in the black hole background. In this case, the metric ( |2.24D becomes the nonrotating BTZ 
black hole: 

£2 



2 2 2 

ds' = -"—JJL de + -dX^ + -A-^ dr^ 



H 



For the usual BTZ black hole, X is written as X = £0 where 
taking X G M. The Hawking temperature ( p.25| ) is, in this case. 



(2.29) 

+ 27r, but here we are 



T = - 

/5 



(2.30) 



In terms of the tortoise coordinate r*, the metric (|2.29|) becomes 



ds' 



H 



-df + drt) + — dX' 



2r 



In 



H 



12.31] 



For the BTZ metric ( |2.29| ), the equation of motion ( |2.28| ) becomes 



H 



X{t,r) 



0. 



(2.32) 



As usual, we proceed by expanding X in modes. Let us set 



X(t,r) 



Ur). 



(2.33) 



^One can show that when the modes on the string arc thermally excited in a black hole background at 
temperature T, this quadratic approximation is valid outside the black hole except for the region within ^/a' 
away from the horizon. 

^In Appendix we will consider the next leading terms (quartic terms) when we estimate the mean free 
path time tmfp- 
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Then the equation of motion ( 2.32| ) can be written as 



,2 , 1 
P 

where we defined dimensionless quantities 



= 0, (2.34) 



p = — , i/= = — . 2.35 

Th Th 27r 

One can see that the hnearly independent solutions to (|2.34|) are given by 
K±) _ 1 P^^^ ( P^\^'''^^ _ _J_p±ii^ +i^r. 



= ^ ^-^ ^ e^^"'^*. (2.36) 



The normahzation in ( p.36 ) was chosen so that, near the horizon, 

fi±) ^ ^±i.r. ^ (2.37) 

and hence the solutions are written naturally in terms of ingoing ("+" sign) and outgoing 
("— " sign) modes. 

2.3 Boundary conditions and cut-ofFs 

Before proceeding with the analysis of the fiuctuations of the scalar field X in the BTZ 
geometry, it is useful to understand the boundary conditions we want to impose on the fields. 
While we are actually interested in the world-sheet theory of the probe string, it is clear that 
we can use the usual AdS/CFT rules to understand the boundary conditions; in the static 
gauge the induced metric on the string world-sheet inherits the geometric characteristics of 
an asymptotically AdS2 spacetime. 

Usually in Lorentzian AdS / CFT one chooses to use normalizable boundary conditions |^ 



for the modes. However, in the present case, that would correspond to a string extending all 
the way to p = oo, which would mean that the mass of the external particle is infinite and 
there would be no Brownian motion. So, instead, we have to impose a UV cut-offQ near the 
boundary to make the mass finite. Specifically, we implement this by means of a Neumann 
boundary condition drX = at the cut-off surface^ 

p = Pc > 1, or r = Tc = VhPc- (2.38) 

The relation between the UV cut-off p = and the mass m of the external particle is 
easily computed from the tension of the string: 

1 r n— _rc-rH [pc - 1) ^Vc 



. dr^/^r = ^:^= ,^ (2-39) 



^°We use the terms "UV" and "IR" with respect to the boundary energy. In this terminology, in the bulk, 
UV means near the boundary and IR means near the horizon. 

^^In the AdS/QCD context, one can think of the cut-off being determined by the location of the flavour 
brane, whose purpose again is to introduce dynamical (and therefore finite mass) quarks into the field theory. 
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Setting 

/.(p) = a[/W(p)+5/(-)(p)], (2.40) 

with constants A and B, we obtain, on implementing the Neumann boundary condition 

dpfuj\p=p^ = 0, 

B = l^ i+Wl ( ^ (2.41) 



1 + iu 1 — ipciy VPc + 1, 

Note that this is a pure phase. This in particular means that, in the near-horizon region 
r ~ rn, we have, because of (|2.37| ), 

X{t, r) = U e-''^' ~ e-*'^(*-''*) + e'^- e-^^(*+"*). (2.42) 

The first term is a mode which is outgoing at the horizon, while the second term is a mode 
reflected at p = pc and falling back into the horizon, with phase shift e*^". The fact that the 
outgoing and ingoing modes have the same amplitude means that the AdS black hole, which 



Hawking radiates, can be in thermal equilibrium at temperature T [p5|] . 

To regulate the theory, we need to introduce another cut-off near the horizon p = 1. 
Specifically, we cut off the geometry by putting an IR cut-off ( "stretched horizon" ) at ps 

p, = l + 2e, e<l. (2.43) 

If we impose a Neumann boundary conditionf^ at ps, we have, just as (|2.41|) , 



B='-^ ^ + + ^ = e-'-(^/^). (2.44) 

l + iu 1 - z(l + 2e)z/ ^ ^ 

If we require ( |2.41|) only, then we determine S as a function of u, but at this point u can take 
any value and is continuous. If e ^ 1, further requiring (|2.44|) effectively makes the possible 
values of u discrete, and the discreteness is given by Au = 1t\ j log(l/e) ^ 1; see Figure ^j. In 
terms of the frequency cu, the discreteness is 

47r2 
/31og(l/e) 

In other words, we have the following density of states: 

All we have achieved by putting the regulator near the horizon is to discretize the contin- 
uum spectrum which naturally occurs when considering horizon dynamics. 



^^One could also take a Dirichlet boundary condition, but in the e ^ limit this would not make a 
difference. 
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arg(B) 




27l/l0g(l/£) 

Figure 2: Boundary conditions at infinity and horizon. First, the UV boundary 
condition ( |2.41|) fixes arg B to lie, say, on the red line; at this point the possible 
values of v are continuous. Further imposing the IR boundary condition (|2.44| ) 
makes the possible values of v discrete (black dots). 



Having done regularization, we can find a normalized basis of modes and start quantizing 
X{t, r) by expanding it in tliose modes. Tliis is standard as in tlie case of a scalar field obeying 
the Klein-Gordon equation; for details we refer the reader to Appendix ^ The upshot of the 
calculation is 



X{t, r) = J2 [(^-^Mt, P) + aiu^{t, pY] , (2.47) 
where the summation is over u discretized according to ( p.45| ). The normalized basis u^^ is 



-ILUt 



(2.48) 



2Puj log(l/e 

where B is given by (p.41| ). The expansion coefficients satisfy the commutation relations 



(2.49) 



2.4 The boundary-bulk dictionary 



Given the behavior of quantum modes on the probe string in the bulk, we can work out the 
dynamics of the endpoint, which corresponds to a test quark in the thermal CFT plasma. 
To understand the precise dictionary we look at the wave-functions of the world-sheet fields 
{X{t,p) in the BTZ geometry) in the two interesting regions: (i) near the black hole horizon 
and (ii) close to the boundary. 

From (|237|) , near the horizon (p ~ 1), the expansion ( p. 471) becomes 



x{t,p^i)^Yl 



a' 13 



2Puj log(l/e) 



iui(t—rf) _|_ giScjg— ia;(t+r») 



) a^ + h.. 



c. 



(2.50) 



We see that the operators are directly related to the amplitude for the outgoing modes, 
g-iaj(t-r*)^ near the horizon. On the other hand, at the UV cut-off p = pc, which we have 
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chosen to be the location of the regularized boundary, ( p.47| ) becomes 
xit)=X{t,pc) = 




2a' (3 



&UJ log(l/e) 



iiy/2 



e-'""' a,, + h.c. 



(2.51) 



This we will interpret as the position of the external particle (test quark) in the boundary 
theory. Here, the operators related to the Fourier coefficients of x{t). 

Using the above relation between (|2.5CI| ) and ( p.51|) , one can predict the correlators for the 
outgoing modes near the horizon, {ciuj-i^alj^ . . .) etc., from the boundary correlators {x{ti) x{t2) . . .) 
in field theory. In particular, if we would be able to make a very precise measurement of Brow- 
nian motion in field theory, we could in principle predict the precise state of the radiation 
that comes out of the black hole. In this way, we can learn about the physics of quantum 
black holes in the bulk from the boundary data. This of course requires us to compute the 
correlation function for the test particle's position in a strongly coupled medium, which is a 
difficult task that we will not undertake here. 

However, at the semiclassical level, we can utilize this dictionary to rather go from the bulk 
to the boundary and learn about the boundary Brownian motion from the bulk data. This 
is possible because, semiclassically, the state of the outgoing modes near the horizon is given 
by the usual Hawking radiation. As argued in [^6|,^, the modes on the string world-sheet 
which impinges on the black hole horizon are thermally excited with a black-body spectrum 
determined by the Hawking temperature. The quickest way to see this is to note that one 
can view our analysis of the fiuctuations of the string world-sheet ( p.27| ) as studying the 
dynamics of massless, free scalars in a two dimensional black hole background.[^ Standard 
quantization of quantum fields in curved spacetime will lead to the modes of the fields 

being thermally excited at the Hawking temperature of this induced world-sheet geometry 
which is the same as that for the BTZ black hole. In particular, it follows that the outgoing 
mode correlators are determined by the thermal density matrix 

Po = 7Tr773?77V' H = 2_^ujala^. (2.52) 



TT{e-PH) ' 



oj>0 



Note that, as we discussed above ( 2.26| ), here we are ignoring the interaction of the string 
with the thermal gas of closed strings in the bulk [vh < r < r^) of the black hole background. 
Namely, we regard the above density matrix ( 2.52| ) as solely due to the interaction with the 
horizon. However, even if we took into account the weak interaction with the thermal gas 
of closed strings, the density matrix would still be given with very good accuracy by ( p.52| ) 
because, at each value of r, the thermal gas is in thermal equilibrium at the local Hawking 
temperature and so is the string. 



^^The induced metric on the string world-sheet clearly has a horizon and is an asymptotically AdS2 space- 
time. 
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As long as we stay in the semi-classical approximation we can use the observations men- 
tioned above to go from the bulk to the boundary and derive the Brownian motion of the 
external particle in the field theory. That is, instead of using the boundary field theory to 
compute the correlation function of the quantum operators a and a\ we can use the fact that 
these correlators are determined by the thermal physics of black holes and utilize them to 
compute the boundary correlation functions. In particular, we propose to use the knowledge 
( p.52|) about the outgoing mode correlators in the bulk, to predict the nature of Brownian 
motion that the external particle on the boundary undergoes. Thus by using the standard 
physics of black holes we will be able to determine the functions 7(t), ^(t) appearing in the 
Langevin equation 



3 Hawking radiation and Brownian motion 

In the previous section, we used the AdS / CFT correspondence to set-up a dictionary translat- 
ing the information about the boundary Brownian particle into corresponding data regarding 
the outgoing modes (the world-sheet oscillators aj and a\ of the fluctuations X^) in the bulk. 
In this section, we explicitly derive the correlation function for the position of the test particle. 
We assume that the outgoing modes are the usual Hawking radiation with the density matrix 
( p.52| ) and derive the result that the endpoint at p = pc ^ 1 indeed undergoes a Brownian 



motion. 

3.1 Brownian motion of the boundary endpoint 

Let us now consider the motion of the endpoint of the string at p = pc ^ 1. We will 
determine the behavior by computing its displacement squared, which corresponds to (|2.3|) . 
In the canonical ensemble specified by the density matrix (|2.52|) , the relevant expectation 
values are given by the Bose-Einstein distribution: 

{al a^i) = Tr (po al a^,) = -J^^- (3-1) 

Using this and (|2.51| ), we compute 

,/N / ,,^/ ^ N> 2a' 3 1 I + 1^^ / 2 cos ut 

.(0)) = ix(t, x(o, pj) = ^ - 



We are using the rescaled frequency u defined in (|2.35|) throughout to avoid clutter. In going 
to the second line, we utilized the density of states determined in (|2.45|) to rewrite the sum 
as an integral. From this, we compute the displacement of the endpoint as: 



s^t) ^ {[xit) - xiO)Y) = ^ / — , coth^ sin2-. (3.3) 
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This has a logarithmic UV divergence. Because this divergence is coming from the zero point 
energy (the e~*'^* term in ( |3.2|) ), which exists even at zero temperature, we simply regularize it 
by normal ordering the a, oscillators (: a^^aj^ : = : a^Oi^ :). When so regularized, the correlator 
fOI) becomes 



:a:txO:=^/ — ^ ^ 7. 3.4 



and the displacement squared (|3l3| ) becomes0'0 



rt) ^ (: - x(0)]^ :) = / - -—^ (3.5) 



We analytically evaluate this integral in Appendix ^. For the present purposes we will only 
record the result for 1, which is all that is relevant for the physics of the boundary field 
theory. We find the following behavior: 

" ^ — t^ (t < tc), 



sUt)-^{^'^: (3.6) 



a! 



t {t:^tc 



So, we observe two regimes, the ballistic and diffusive regimes, exactly as for the standard 



Brownian motion ( |2.3| ). The crossover time tc is given by 



a' m 



tc-/3pc--^,. (3.7) 



In the ballistic regime, t -C the coefficient of t in ( |3.6|) is exactly the same as (|27 
determined by the equipartition of energy x ~ ^JT jm. In fact, one can say much more; in 
Appendix 0, we show that, if pc ^ 1, the probability distribution /(p) for the "momentum" 
■p = mx oi the endpoint is exactly equal to the Maxwell-Boltzmann distribution for non- 
relativistic particles, 

/(p)oce-^^^ Ev = ^- (3-8) 



^'^Note that regularization by normal-ordering docs not preserve the KMS relations except in the classical 
limit. 

-'^^ Another way to regularize the correlator is to use the canonical correlator introduced in ([4.7]). Using 



(4.12), one can derive 

sl{t) ^ l\x{t) - x(0)]; {x{t) - x(0)]) = ^ / ^ sin^ 



2a'/32 r°° duj 1 + . ^ tot 



(K|//?)-(l-p-2)(l-e-2-l*l/^''0' 



This is finite and has exactly the same short- and long-time behaviors as in ( |3.6| ). Note, however, that 



the divergence ( |3.3[ ) is related to the well-known fact that the fluctuation of the position of a string always 
diverges |3S|]. 
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In the diffusive regime, t ^ tc, we find a diffusion constant (one half of the coefficient of 



a' 



^AdS3 = ■ (3.9) 

A priori this looks counterintuitive because it is inversely proportional to temperature T and 
implies that the random walk becomes more vigorous at lower temperature. However, this 
is consistent with the known results for test quarks moving in the thermal A/" = 4 super 
Yang-Mills plasma [^,|2l|,|2^,|24]. For example, refs. [pOj^ considered a heavy particle on 
the boundary moving at a constant speed v under the influence of an external force. One can 
compute the friction acting on it from the bulk using the Nambu-Goto action, where a string 
is moving at velocity v, trailing along the boundary. It is easy to generalize their computation 
to AdSrf with general d, the result bein^ 



^ ~ " (rf^^Tj^ (r^"^;^jv(5^' p-^fT^- ^^-^^^ 

In the non-relativistic limit, f <^ 1, this means that the friction constant is 

If we use the Sutherland-Einstein relation ( p. 61) 1*^ we obtain the diffusion constant 

which agrees with ( |3.9| ) for d = ?>. 

One can give an intuitive explanation for the reason why the diffusion constant is inversely 
proportional to T from the boundary viewpoint of Brownian motion. The random walk 
behavior of Brownian motion is due to frequent collisions of the Brownian particle with 
the fluid particles. In particular, after n steps (collisions), the distance s that a random- 
walk particle covers scales as y^Lmfp, where the mean free path Lmfp is the typical length 
traveled between the collisions, i.e., it provides a scale for the system. For the thermal system 
under consideration we have, L^fp ~ 1/T, because this is the only scale available in a CFT 
at temperature T.0 After n collisions, the time elapsed is given by t ~ n/T., since the 
time between collisions is also given by L^fp ~ 1/T. So, putting things together, we have 
s ~ y/tT ■ 1/T = ^Jt/T, namely, ~ t/T which is exactly what we infer from ( p.l2| ). 



Although this is a straightforward generahzation of [|0|,^ and the general formalism has been laid out 
in |p3| , it seems to us that this result has not appeared explicitly in the literature. 
^^Note that, as explained around (2J), the relation ( |2.6| ) does not depend on d. 

^®The precise value of the Lmfp depends on the strength of the field theory couphng, but the temperature 
dependence follows via dimensional analysis. In fact, in Appendix we estimate the mean free path time 
to be trat-p l/(V^r), where A ^ £*/q;'^. If the plasma constituents are moving at the speed of hght, this 
means that i,nfp ~ l/{-/\T). With this value of L^fp, we can even recover the A dependence of ( ^.12 ). 
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From the bulk point of view, on the other hand, one can give a physical explanation for 



D ~ 1/T as follows. Near the horizon (p ~ 1), the Nambu-Goto action (|2.27|) becomes 

^NG ^ I dtdr^ [{dtXf - [drXf] . (3.13) 

This is the same as the action for a string in flat space, with a' replaced by a^^j = a' / (4-7r^^^T^). 
This means that the size of the fluctuations in X"^ is proportional to a^fj ~ T~^. The Boltz- 
mann factor of Hawking radiation gives an additional factor of l/{e^^ — 1) which scales as 
T at low frequency. Altogether, near the horizon, the fluctuations scale with temperature as 
~ T^^. When a fluctuation propagates to p = pc, a greybody factor damps the fluctu- 
ation. However, as one can see from ( p.5lD , the damping is 0{1) for very small frequency. 



This leads to = X{p = pcY being ~ T~^. The reason why very low frequency modes can 
reach p = undamped is that X is an isometry direction and very low frequency X modes 
can propagate at almost no cost in energy. 

A natural question to ask is what happens to this T"^ scaling as T ^ 0, as we expect that 
the endpoint should not fluctuate at T = 0. This can be understood by realizing that for a 
mode to propagate to p = pc undamped, it should be the lightest mode in the problem. In 
particular, only modes whose frequencies are lower than the thermal scale (which goes to zero 
as T — i> 0) can propagate without damping. Translating to real time dynamics, this means 
that one needs to wait until t ~ to see the diffusive regime; but since oc — oo as 
T — > we never enter that regime and the motion is always ballistic as expected. 

Thus, we have demonstrated that the endpoint of the string at p = pc indeed behaves like 
a Brownian particle; it shows ballistic and diffusive regimes, just as for the usual Brownian 
motion. We would now like to understand the Langevin equation from the bulk perspective. 
As we will see below, the Langevin equation governing this Brownian motion turns out to be 
not of the simplest type ( |2.1| ) and ( p.2| ), but rather the generalized one ( |2.8| ) indicating that 
the precise nature of the random kick encountered by the Brownian particle depends on the 
past history of its trajectory. 

3.2 Forced motion & the holographic Langevin equation 



As we discussed in subsection |2.1| , the generalized Langevin equation ( p.8| ) has two functional 
parameters: the memory kernel 7(t) and the auto-correlation function K,{t), related to the 
dissipative and stochastic components, respectively. We would like to determine these func- 
tions from the holographic viewpoint for the probe string in the black hole background. In 
order to do so, we will first determine 7(t), or equivalently p{uj). Once we know p{uj), we can 
compute k{uj) by using equation ( p.l9| ) and the {xx) correlator (|3.2| ) (or (|3.4|) ). 

We first turn to the determination of p{uj). Consider applying an external force on the 
Brownian particle as in (|2.14|) ; from the response to this force we can read off p{uj) using 
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( p. 15 ). So the natural question is what external force is to be applied to the string endpoint. 



As in the AdS/QCD set-ups, we can realize such forced motion by placing a "flavor" D-brane 
at the UV cut-off p = and by turning on world- volume electric field on it. Since the 
endpoint of the string is charged, this will exert the desired force on the Brownian particle. 

So, let us consider the Nambu-Goto action ( 2.27| ) in the general metric ( p. 26 ), and add to 
it the following boundary term 

5bdy = ^A(a:,X), (3.14) 

which corresponds to turning on world-volume field on the flavor D-brane (which is placed at 
the UV cut-off pc). Here, A{x, X) is a 1-form defined on the flavor D-brane world- volume. We 
again work in the gauge where the world-sheet coordinates are identified with the spacetime 
coordinates = t,r. We have = X\x) and 

Sbdy = ^[Atix,X)+Ai{x,X)X']dt, (3.15) 

where t is taken to be the coordinate along the boundary as before (or equivalently, the 
boundary is at r = const.). The equation of motion one obtains for the total action Sj^jq -|- S'bdy 
at the boundary is 

n'^Gjjd^X' - 2W {Fjt + FjjdtX') = 0, (3.16) 

where g^i, is the induced metric on the boundary, is the outward-pointing unit normal to 
the boundary, and Fjt = djAt — dtAj, Fjj = djAj — djAj. 

Returning to the simple setting of the BTZ geometry (|2.29|) , the equation of motion for 
the string in the presence of this additional gauge field is 

p^{p^-l)d,X = ^^Fj,, at p = p,. (3.17) 

For the world-volume field Fxt we choose an oscillating electric field with frequency u: 

Fxt = E = E,e-'^\ (3.18) 



motivated by ( p.l4| ). We now want to compute how the string, in particular its endpoint 



X[t,Pc) = x{t), moves under the influence of this external force in order to compute the 
admittance /i(cij). 

As before, the solution to the bulk equation of motion can be written as a linear combi- 
nation of the modes /^^''(p): 

Xit,p) = [A'fi-'\p)+B'fi-\p)] e--*. (3.19) 
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To determine the coefficients A', B', we need to impose a boundary condition at the horizon, 
in addition to the boundary condition ( p.lTp at p = pc. Effectively all that the external field 
has done was to modify the Neumann boundary condition, which we imposed earlier on the 
cut-off surface, to a mixed boundary condition. 

In the semiclassical approximation, the boundary condition near the horizon is such that 
outgoing modes are always thermally excited because of Hawking radiation, while the ingoing 
modes can be arbitrary.]^ From (|2.37|) , the coefficients A' and B' correspond to outgoing and 
ingoing modes respectively. Therefore, the boundary condition at the horizon is that A' is 
thermally excited. However, because the radiation is random, the phase of A' takes random 
values and, on average. A' vanishes: {A') = 0. Recall that it is such averaged quantities that 
we are interested in; the admittance fi{uj) is obtained by suitably averaging over the ensemble, 

cf. (in. 

Requiring the boundary condition (|3.17f) , with electric field (|3.18| ) at p = pc and the 
condition that (A') = 0, we determine the average values of A', B' to be 

(A') = 0. m = ''"f n"'\ { ^) Bo ■ (3.20) 

r}j z/(l - ipcv) VPc + 1/ 

From this we infer that the average value of X at the UV cut-off p = pc is 

(.W)HA-«,P.))^^J^^^;„e-, (3.21) 
which in turn implies that the average value of the momentum p = mx is 

Comparing this with ( p.l5| ), we obtain the admittance 

, , 1 a' B"^ m 1 — iu / Pc 
n(uj) = — = ^1- 'J^. 3.23 

A simple check on the consistency of these computations is to compute the energy flow 
along the string falling into the horizon, which must be equal to the work done by the external 
force. For the theory ( p.27|) , the stress-energy tensor is 

= {a''' ^7"') Gij d^X'd.X'. (3.24) 

Zn a \ I J 

Because we are working in static gauge this world-sheet stress-energy tensor measures the 



spacetime energy. In the case of the BTZ spacetime (|2.29| ), the flow of energy along the r 
direction is 

= r''{r''-r\) j^^r^ ^ 9^X1 = RefSpXS^X]. (3.25) 



^^If one ignores the thermal excitations of the outgoing modes and sets them to zero, this boundary condition 
becomes the so-called purely ingoing boundary condition. 
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Here, we replaced d^XdxX Keld^XdxX] so as to work directly with complex fields. Con- 
sider the solution for X{t, p) as in (|3l9|) with the coefficients A', B' given by the average value 
( p.20|) (we ignore thermal fluctuations in replacing the amplitudes by their average). Then, 
(PD evaluates to0 



271 a' fE^ 1 + 1/2 



V^T[ = 2 - ,,00 • (3-26) 

On the other hand, the work done per unit time (namely, power) by the electric field E acting 
on the endpoint at X(t,pc) is 

Re[EdtXit,p,)] , (3.27) 



where E is given by ( |3.18| ). For X{t,p) in ( |3.19| ), it is easy to check that this equals ( |3.26| ). 



Hence indeed as expected, the work done by the external force is transmitted down the string 
into the black hole horizon and the energy is thus dissipated away. 

3.3 The holographic auto-correlation function and time scales 

We now turn to the computation of the random force correlator, k{uj). From the (xx) 
correlator (|3^), we can compute the (pp) correlator as 

, , s , s , r^o , / N / N > a' 0"^ m? /"°° , 1 + 1/2 ^ 

{■.pit)piO):) = -m'dU:xit)xiO):) = IT^ 

(3.28) 



7r£2 27r 1 +p2;y2 g/3|<^| _ ^ " 

To obtain the power spectrum defined in ( |2.16| ) for the momentum p we Fourier transform in 
time t to obtain 

Here, the superscript "n" is for remembering that this power spectrum was computed using the 
normal ordered correlator {:pp:). Then we can exploit the relation ( |2.19| ) between the power 



spectrum for the auto-correlation function and the momentum spectrum and the previously 
derived expression for p{uj), ( |3.23| ), to obtain the power spectrum for the random force R, 
which is nothing but the random force correlator k'^{uj): 

.-H,A r(,A 47rf2 l + z/2 



^"^The plus sign is because t is a lower index. If we raise t, this will have a minus sign, indicating a flow of 
energy toward the direction of the horizon (smaller r). 
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Next, let us compute the physical time scales treiax and tcou- First, from (|3.23|) , one can 
compute the relaxation time treiax defined in ( |2.20 ) as: 

irelax = = 0) ~ — . (3.31) 

To compute the collision duration time tcou, we first need the real time auto-correlation 
function for the random force (RR): 

n-{t) = (: R{t)R{0) :) = j^^ ^ I^u^) e^*. (3.32) 

By using the explicit form (|3.30|) , we obtain 

'^"W = 3^ [plh,it,P) - {pi - 1) /i2(t,Apc)] , (3.33) 
where we defined the functions 

/oo ^-itx/p roc 11 ^-itx//3 

eN_i ' ^2(t, P, Pc) = J Jx ^ _ (3-34) 

For Pc ^ 1, hi and /i2 are almost equal; if x <^ pc, we can approximate 1 + {j^Y the 
integrand of h2 by 1 while, if x > pc ^ 1, the integrand is almost vanishing because of the 
Bose-Einstein like factor l/(e'^' — 1). Therefore, the RR correlator evaluates to 



2r , , 2£' 



t J sinh'(7rt//5). 



(3.35) 



This function has a support of width of order (3 around t = 0. Therefore, using ( p.22| ) we 
obtain the collision duration time 

tcow ~ /3 = ^, (3.36) 

The T dependence is as it should be from dimensional analysis in a CFT at temperature T, 
but the fact that this is independent of the 't Hooft coupling A is not trivial. 
The ratio of the two time scales is given by 

^ ^ ^ ^ ^ (3 37) 

where we related a' / P to the boundary 't Hooft coupling]^ by using the relation ~ A 

[^, In the weak or moderate coupling regime, A < 1, we can make this ratio large 
by considering a Brownian particle with m ^ T and obtain the standard Brownian motion 
as explained below (|2.23|) ; the Brownian particle becomes thermalized only after numerous 



^^Strictly speaking, this is only the 't Hooft coupling A in the standard AdSs case, but we will use the same 
terminology to denote i'^/a''^ for other values of d as well. 
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collisions with fluid particles. In the strong coupling regime, A ^ 1, however, this is not 
the case and, in order to have the standard picture, we have to consider a much heavier 
Brownian particle with mass m ^ \/~\T , which is always possible. On the other hand, if 
T <^ m <ti \fXT, the situation is totally different. Although the effect of a collision with a 
single fluid particle (with energy ~ T) is small, because the Brownian particle interacts with 
many fluid particles at the same time, it can become thermalized in a time much shorter than 
the time it takes for a single process of collision. To make this claim more quantitative, we 
estimate in Appendix the average time tmfp between collisions. The contribution that a 
single collision makes to the random force R{t) has width tcou- R{t) consists of many such 
contributions, with the typical distance in time between two collisions being tmfp- Determining 
tmfp is not entirely straightforward, as it requires us to analyze the four-point correlation 
function of the random force, and we only flnd a non-trivial answer once we take the fourth 
order correction to the Nambu-Goto action into account. As a result, t^aip is suppressed by a 
factor of l/\fX compared to tcoih the flnal result being 



At weak coupling A <^ 1, we have tmfp ^ ^cou and the standard kinetic theory picture, where 
the Brownian particle is occasionally hit by a fluid particle, is valid. On the other hand, at 
strong coupling A ^ 1, we have tmfp ^ ^cou, namely, many collisions occur within the time 
scale for a single collision process to take place.^ This supports the picture above that the 
Brownian particle interacts with many fluid particles at the same time. 

4 Fluctuation-dissipation theorem 

Thus far we have seen how the string probe in the bulk geometry holographically captures the 
Brownian motion of an external test particle introduced in the boundary CFT plasma. As we 
have seen explicitly, one can derive the Langevin equation for the string endpoint by tracing 
back the information about the part of the string that is touching the black hole and hence 
gets thermally excited due to the outgoing Hawking quanta. One of the hallmarks of non- 
equilibrium statistical mechanics is the fluctuation-dissipation theorem which relates 
the observables in the system perturbed inflnitesimally away from equilibrium to equilibrium 
quantities. We now turn to show that not only are the results we derived in section ^ consistent 
with the fluctuation-dissipation theorem, but that we can in fact obtain this result directly 
from the gravity side. 

^^Perhaps the term "mean free path time" is not an appropriate one in this regime where a second colhsion 
takes place before the first one ends, and thus the particle is never freely moving. However, there being no 
other choice, we will continue to use this term in the strongly coupled regime. 
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4.1 Linear response theory 

We begin our discussion of the fluctuation-dissipation theorem with a hghtning review of 
hnear response theory ^ . 



Consider a system whose unperturbed Hamiltonian is given by H . Assume that, in the 
inflnite past t = — oo, the system was in an equihbrium state with the density matrix 

P' = (4-1) 

Now perturb the system by adding an external force K{t) conjugate to a quantity A. The 
total Hamiltonian is 

i/tot = H + H,,,{t) = H -AK{t). (4.2) 

Under this perturbation, the change in another quantity B is given, to the flrst order in the 
perturbation -ffcxt, by the so-called Kubo formula: 

AB{t)= I dt'K{t')<f)BAit-t'), <PBAit)^~t{[A{0),B{t)]), (4.3) 

where we deflned (O) = tr(peC) and 0{t) = e^^^Oe~^^*. The function (pBAi^) is called the 
response function. 

If we consider a periodic force with frequency u, 

Kit) = Ko e"*^*, (4.4) 

then ( [4.3|) gives the following change in B: 

AB{t)=i^BA{oo)Koe-'''\ (4.5) 

where the admittance ^iba^^) is given by 



1 



oo 



fiBAiuj)= dt<PBA{t)e'^' = - dt{[A{0),B{t)])e''^' = /3 dt {AiO); B{t))e'-\ (4.6) 
Jo ^ Jo Jo 

and with the canonical correlator (X; Y) deflned by 

(X; y) = -p[ d\ {e'^Xe-'^'Y) = \ (^(-^A) Y) , (4.7) 

which satisfles the following properties 

(X(0); Y{t)) = {Y{ty, X(0)) = (r(0); X{-t)). (4.8) 

The relation (|4.6|) is called the fluctuation-dissipation theorem, because the right hand side 
is the fluctuation (correlator) in the equilibrium state pe, while the left hand side yields the 
admittance which is related to the dissipation (friction). 
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In the case of Brownian motion, we can take A = x and i^ext = —xK{t), where K{t) is 

identified with the external force appearing in the Langevin equation (|2.^ ). Then, for B = p, 
we obtain the admittance 

/xM = ^ r dt{p{Q)-p{t))e^-\ (4.9) 
Jo 

Due to the relations (|4.8|), this implies 

2Re/iM = - r dt{p{Oy,p{t))e'-' = ^I^iu), (4.10) 

where Ip{uj) is the power spectrum for p defined using the canonical correlator. From this, 
using the relation ( |2.19| ), one can derive a more direct relation between the friction and random 
force as 

2Re-f{uj) = ^lUu) = ^K'^iu), (4.11) 
m m 

which is sometimes called the second fiuctuation-dissipation theorem, in contrast with 
or (|4.10|) which is sometimes called the first fiuctuation-dissipation theorem [^,^ . 



4.2 Explicit check of fluctuation-dissipation theorem 

The fiuctuation-dissipation relations ( [4.9| ), ( [4.10| ), and ( [4.11| ) for Brownian motion were de- 
rived from the field theory viewpoint and are not immediately obvious from the bulk view- 
point. Here, let us explicitly check that they indeed hold using the explicit results obtained 
from the bulk in section |^. 

Similarly to ( p. 21) or ( p.4|) , we can compute the canonical correlator for x as 



W0);.W)^^r^!^A^-" (4,12) 



27r 1 + pIu^ 

Because p = mx, this implies the following canonical correlator for p: 



m-Mt))-"^r^^e--'- (4,13) 



■kP 27r 1 + p2iy2 

This means that the power spectrum for p is 



On the other hand, from (|3.23| ), we immediately obtain 



2Re/xH = 1^22 - 4.15 
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By comparing (|4.14 ) and ( |4.15|) , we see that the fluctuation-dissipation theorem of the form 



( 4.10 ) indeed holds. This also implies the second fluctuation-dissipation theorem ( 4.11|) 



These relations can be regarded as providing evidence that the motion of the string end- 
point in the bulk can be described by a generalized Langevin equation. Note, in particular, 
that the way we derived the fluctuation (correlator) and the way we derived dissipation (ad- 
mittance) were very different; for the former, we assumed thermal Hawking radiation near 
the horizon and measured the position of the string endpoint at the UV cut-off while, for the 
latter, we considered forced motion imposing a boundary condition at the horizon which was 
essentially the purely ingoing boundary condition. In the next subsection we describe how 
these two quantities are related directly from the bulk point of view. However, it would be 
desirable to have a more intuitive physical understanding of why this should be the case. 

4.3 Bulk proof of fluctuation-dissipation theorem 



In subsection |4.2| , we demonstrated that the fluctuation-dissipation relations holds for the 
special case of string probes in the BTZ spacetime by an explicit bulk computation. We 
now prove that the fluctuation-dissipation relations hold more generally, again from the bulk 
viewpoint. 



Consider a string probe in the d- dimensional metric (^.2(j|). We would like to turn on an 



electric field = Ej{t) on the flavor D-brane aX r = and consider the resulting position 
x\t) = X^{t,rc) of the string endpoint in response to it. If we take At = Ejit) X^, Aj = 0, 
then the boundary action, ( p. 151) , can be written as 

^bdy = jdtEi{t)X^ = J dtdrSir -r^)Ej{t)X^. (4.16) 

This can be regarded a source term for the field X^; upon inclusion of this term, the equation 
of motion ( p.28| ) is changed to 

[Gjj{x) d,X\x)] =-^6{r- r,) Ej{t). (4.17) 

As is standard, we can solve this by using the retarded propagator 

D^J,{t, r\t\ r') = e{t - t'){[X\u r).X-\t', r')]), (4.18) 



where X\t,r) satisfies the equation of motion ( |2.28| ) (or equivalently (|4.17| ) with the right 



hand set to zero) and can be expanded in modes as in ( |2.47| ). Namely, 

X\t,r) = Y,K{t,r)a^ + ui{t,rraU, 

oj>o (4.19) 
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where {ulj{t,r)} is a normalized basis of solutions to ( p.28 ). Note that, with ([4.19|) , the 
commutator appearing in ( [4.18 ) is actually a c-number and -D^t is independent of the state 
with respect to which we take the expectation value. D^^^ can be shown to satisfy 

[Gu{x) d.Di^it, r\t\ r')] = ^ ^ - t') Sir - r'). (4.20) 

Therefore, the solution to ([4.17|) can be written, using -D^i, as 

X\t,r) = i / dt'D'J^{t,r\t',r,)Ej{t') 



ret V 

-oo 
oo 



(4.21) 



dt"{[X'{t,r),X'{t-t",r,)])Ej{t-t"). 



By setting r = Tc, we obtain the position x\t) = X\t, Vc) of the string endpoint in response 
to the external force Ej{t) as: 

POO 

x\t) = % \ dt" {[x\t),x^it - t")])Ej{t - t"). (4.22) 



If the fluctuation-dissipation theorem is to hold, this must be equal to the x^{t) obtained 
by using the Kubo formula (|4.3|) , identifying x\t) with the position of the Brownian particle 
in the boundary. From the action ( [4.16|) , one reads off the external force appearing in 
to be AK = Ejx"^ . Then, by applying the Kubo formula ([4.3|) for B = x^ , 



x\t) = -i I dt' Ej{t'){\x\Q),x\t-t')]) 

(4.23) 



oo 
oo 



dt" {[x\t"),x\mEj{t-t"), 

where t" = t — t'. Because the system is stationary, the expectation value is invariant under 
shift of time: {[x\t"), x-^ (0)]) = ([x^(t), x^(t - t")]). Therefore, the bulk response ( ^2% is 
the same as the expression (|4.23|) computed from the boundary Kubo formula. This implies 
that the fluctuation-dissipation relations ( |4.9|) , ([4.10|) , (|4.11|) indeed hold. 



5 General dimensions 

Thus far we have considered the case of d = 3 dimensional asymptotically AdS spacetimes 
only, which had the advantage that the wave equation for the modes of the string was exactly 
solvable. For general d, this is no longer possible and we have to use approximate methods. 
In this section, we employ the low frequency approximation u <^ T and briefly summarize 
how some of the results of the previous sections get modified in the case of asymptotically 
AdSrf spacetimes with general d. 
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The starting point is the metric (|2.24 ), with Hawking temperature given in ( 2.25 ). The 
tortoise coordinate r* is defined via 



1)2 



dr^ = dr. (5.1) 

n[r) 

If we define rj = exp[27rz/ (d — 1)], an exphcit expression for the tortoise coordinate is 



k=0 



The term with k = shows that near the horizon, r^, ~ ^^J^^ log(^ — l), and from (|57 
we also see that the behavior near infinity is 

e 

~ , r ^ oo. (5.3) 

r 

The generahzation of (p.32|) to arbitrary d reads 

- dlX + ^ dr[r^ h{r) drX] = 0. (5.4) 

As before, we will exploit the translational invariance along t to decompose modes in plane 
waves; for convenience consider solutions of the form 

X(t,r) = e-^'^*r"i<l><^(r) , (5.5) 

from which it follows that the functions ^ui{f) satisfy the following equation 



dr 



2 +u^-V{r) 



^u{r) = , (5.6) 



with 



V{r) = j^r^ h{r) [2 h{r) + r h'{r)] . (5.7) 



The wave equation (|5.6| ) can be thought of as a time-independent Schrodinger equation for a 
particle moving in potential V{r). 

As in section ^ we want to exploit the semiclassical physics of Hawking radiation to learn 
about the behaviour of the string endpoint on the boundary. Once again it is worth noting 
that the dynamics of the scalar field X(t,r) is similar to a minimally coupled scalar field 
propagating in an asymptotically AdS2 spacetime with an event horizon. We would like to 
compute the admittance for the Langevin equation in this case. 

In order to redo the computation that led to (|3.23|) we need to find the solution of the wave 



equation (|5.6| ) which is purely ingoing at the horizon r = rn- Let us denote this particular 
solution of the wave equation by X~{r). It is not possible to obtain this for general frequencies 
(jU and hence we employ a low frequency approximation u <^T and use the so-called matching 
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technique. Here, we only write down the final result of the computations, relegating the details 
to Appendix |^. The solution that is purely ingoing at the horizon behaves near infinity as 



Xj(p) = C+X+(p) + C"X^(p), (5.8) 

where 

C^ = \ (l±^ + ^6^), X±(p)= (^iT^^e^-/''. (5.9) 

Here, 6 is a constant independent of z/, whose precise value is not relevant for our purpose. 
Also, as before, we defined dimensionless quantities 

p = — y = — • 5.10 



In terms of these, the low frequency condition uj <^T reads u <^1. Actually, this result (|5.8| ) 
is valid only to leading order in the v expansion. 



By carefully redoing the calculation in subsection p.2| , one can show that there is the 



following relation between the ingoing mode and the admittance ^{oj): 



Using the explicit expression of (Eqs. (|0|), (|5.8| )), the final result is 

_ i{d-lfa'ml3^ [{l+tbu)up,-t] v cos(^) +z \p,- iv^ {\ ^ ihv)\ sin(^) 

(5.12) 

As mentioned above, this result is valid only to leading order in the expansion in v. For small 
z/, the right hand side of ( |5.12| ) behaves as 



, , ((i — 1)^ m , 



For c? = 3, this agrees with (|3.23| ) for z/ ^ 0. Furthermore, this agrees with the drag force 
result ( p.ll| ) for general d, because 70 = p[0]~^. 



Just as we did for the ci = 3 case in subsection 3J^, we could also compute k{uj) for 
general d in the low frequency approximation. However, this is not necessary, because we 
can directly obtain from the fiuctuation-dissipation theorem ( 4.11| ), whose validity we 



already demonstrated for all values of d in subsection |4.3| . 

6 Stretched horizon and Brownian motion 

The main philosophy of the membrane paradigm is that, as far as an observer stayin; 



outside a black hole horizon is concerned, physics can be effectively described by assuming 
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that the objects outside the horizon are interacting with an imaginary membrane, which is 
endowed with physical properties, such as temperature and resistance, and is sitting just 
outside the mathematical horizon. In section ^, we assumed that the Brownian motion of the 
UV endpoint of a string was caused by the boundary condition we impose at the horizon — all 
ingoing modes are falling in without being reflected, while the outgoing modes are always 
thermally populated. A curious question then is whether this boundary condition can be 
reproduced, in the spirit of the membrane paradigm, by postulating some interaction of 
the string with a membrane at the stretched horizon just outside the actual horizon. The 
interaction necessarily assumes a stochastic character, so it is natural to expect it to be 
described by a sort of Langevin equation. For a schematic explanation, see Figure |^. It must 
be noted that the physics of the stretched horizon has been discussed in the AdS / CFT context 
previously in [^^,|^,^ and more recently in [^| where there is a nice discussion regarding 
the dynamics of the stretched horizon and the universality of hydrodynamic coefficients. We 
will now turn to a derivation of the properties of the stretched horizon in section |6.1| and 
then proceed to ask whether we can learn anything about the microscopic structure of the 



stretched horizon in section 6.2 



6.1 Langevin equation on stretched horizon 



Brownian 
motion 



UV cutoff 






string 



■p=p, 
p 



•X 



stretched 
horizon^ 



mathematical^ 
horizon 



force acting 
on IR endpoint 

p=\+2z 

-p=\ 



Figure 3: A membrane-paradigm like picture of the Brownian motion. There are 
friction and random force acting on the IR endpoint of the string on the stretched 
horizon, effectively giving the boundary condition. 



Let us consider placing an imaginary "IR brane" near the horizon atps = l + 2e,e<^l 
and assume that the string ends on it.0 If we assume that a force is acting on the endpoint, 

^^In reality a string that dips into the black hole will continue merrily past the horizon without any trouble; 
the quickest way to see this of course is to pass to coordinates that are regular on the horizon such as ingoing 
Eddington-Finkelstein or Kruskal coordinates. Here we are interested in mimicking the boundary conditions 
of the black hole and hence will postulate there to be an imaginary boundary in the IR at p = p^. 
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the equation of motion for the endpoint is, just as in (p.l7|) , given by 



-f^W. = ff. (6.1) 

where is the force along the X direction measured with respect to the time t. Note that 
there is no term hke mX on the left hand side, because the endpoint has zero mass, having 
zero length. We assume that the force , just as in the usual Langevin equation (^^ ), has 
frictional and stochastic components: 

F^{t) = -f dt'^,{t-t')dtX{t',ps) + Rs{t), (6.2) 

J — oo 

{Rs{t)) = 0, {Rs{t)R,{t')) = ^s{t - 1'), (6.3) 

where we allow the friction to depend on the past history through a memory kernel 7^. We 
would like to choose 7^ and Kg appropriately to reproduce the correct boundary condition 
described above. 

Near the horizon the fluctuation of the string is given by (|2.50| ) 



where on is discretized with Auj given in a^A^ . aL+^ and annihilation operators 

for outgoing and ingoing modes, respectively. Depending on the boundary condition one 
imposes at the UV cut off, and ai~^ get related to each other (for example, in the 
case of the Neumann boundary condition we imposed in subsection |2.3| , they are related as 
However, because we are considering a Langevin equation which holds 
independent of such relations, we regard and independent variables. 

Plugging ( |6.4|) in and going to the frequency space, we can write the equation of motion 
rri) as 



a' Pu 



2& log(l/e) 



Rs{u;) 
(6.5) 



for u; > 0. Here, 7s [c*;] is the Fourier-Laplace transform of 7s (t) similar to (|2.12| ) while Rs{uj) 
is the Fourier transform of Rs{t) as in ( p.ll|) . In order to realize the boundary condition that 
all ingoing modes fall in without reflection, we should set the coefficient of aL ^ to zero (since 
we want to be able to set the ingoing amplitude ai""* to any value). This gives 

^■M = 2^^ = ^7]P ^ 7.W = ^«W. (6.6) 
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Substituting this back into (|6.5| ), we obtain the relation between the random force and the 
outgoing mode coefficients 



If this random force is to reahze the thermal nature of the outgoing modes, {a^^'^a'^^) = 
[e^^ — then from ( |6.7|) we obtain 

{RM'RM) = (6.8) 

for/3tu<l. (6.9) 



a' (3^ ' 

Here 'D{ijj) is the density of states defined in (p.46|) . This means that the correlator for the 
random force is 

dujViiu) (i?,,(^)t R,{uj)) e'-(*-*') ^^5{t- t'). (6.10) 

-oo ^ P 



The delta function behavior is due to the approximation we made in (|6.9|) ; the actual Kgit — t') 
is nonvanishing for |t — t'| ^ /5, as one can see if one uses the original exact expression (|6.7|), 



In summary, the boundary condition near the horizon can be effectively realized by a 
Langevin equation for the string endpoint X{t,ps) at the stretched horizon given by0 



^Jf. 5pX = --^ d,X + RM, {RMRsit')) - ^ ^(^ - t'). (6.11) 

The two terms on the right hand side of the ffist equation are, respectively, i) friction which 
precisely cancels the ingoing waves, and ii) random force which is responsible for the outgoing 
modes being thermally excited at the Hawking temperature. 

Given the auto-correlation function for the random force Rs{t) acting on the string end- 
point at the stretched horizon, we can exploit the Sutherland-Einstein relation ( p.4|) to com- 
pute the diffusion constant on the stretched horizon. We find 



2T2 _ a' 
K,(0) ~ 2TrPT 



A. \Jc 

^IdSa = -777^ = 7737^ (6-12) 



which is the same as the diffusion constant for the string endpoint undergoing Brownian 
motion in the boundary (|3.9|) . In deriving (|6.12| ) we had to assume that the dynamics of the 
string endpoint on the stretched horizon obeys the Sutherland-Einstein relation derived for a 
point particle. In other words, we assumed that a point particle fixed on the stretched horizon 



*Note that these relations are operator relations whose full structure has been given in (bj) and (WM 
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external 
string 




Figure 4: A possible microscopic picture of a black hole, where the horizon is 
covered by a stringy "cloud" made of strings and branes. An external fundamental 
string ending on a horizon is dissolved into the cloud and incessantly kicked around 
by the cloud. 



will experience the same friction and random force as the ones appearing on the right hand 
side of ( |6.11| ), and thus will random walk with the diffusion constant (|6.12| ). In fact, we will 



now argue that this is not quite unexpected from the viewpoint of the membrane paradigm. 

In the context of the membrane paradigm, it is conventional to ascribe transport properties 
to the stretched horizon. In fact, it is well known that the shear viscosity of the black hole 
membrane saturates the famous bound derived in the boundary field theory, rj/s = l/47r, 
cf., [^,^.0 In the hydrodynamic regime of the AdS/CFT correspondence, Ref. |^6| argued 



that one can derive the universality of this ratio using the membrane paradigm, i.e., the 
physics of the stretched horizon similar to the discussion given above. We have here focussed 
on the stochastic Langevin process and derived the features of the membrane that reproduce 
the physics of strings impinging on the black hole. Again we see that the diffusion coefficient 
of heavy quarks in the boundary ( p.9|) agrees with that derived for the stretched horizon 

(in. 

6.2 Granular structure on the stretched horizon 



In P9|,|50[], Susskind and collaborators put forward a provocative conjecture that a black hole 



is made of a fundamental string covering the entire horizon. Although this picture must 



be somewhat modified |51| since we now know that branes are essential ingredients of string 
theory, it is still an attractive idea that, in the near horizon region where the local temperature 
becomes string scale, a stringy "soup" or "cloud" of strings and branes is fioating around, 
covering the entire horizon. 

If this picture is true, one naturally expects that it is this stringy cloud that is exerting 
frictional and stochastic forces on the IR endpoint of the fundamental string as described by 



25c 



See also [48| for another membrane paradigm inspired perspective on the ratio 77/s. 
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( 16.11 ); see Figure^. Can we learn anything about this stringy cloud? The stretched horizon 



is located a distance ~ = va' away from the mathematical horizon It is occupied by 
a string of length L ~ 5*/^, with S the entropy of the black hole. If we associate one degree 
of freedom to each string segment of length Is, the number of degrees of freedom equals the 
entropy, and one can try to think of these degrees of freedom in terms of free quasi-particles. 
The average separation between the quasi-particles is equal to 

AX~— (6.13) 

with /p(rf) the (i- dimensional Planck length. If the quasi-particles move with the speed of light 
and a represents the probability that quasi-particles will interact with the endpoint of the 
string, then we expect a mean free path time of the order of 

AX 

tmfp . (6.14) 

a 

Supposing for a moment we assume that this is the same as the mean free path time on 
the boundary which, in Appendix 0, we argued to be given by t^{p ~ 1/(T-\/A). Combined 
with ( |6.14| ) and (|6.13|) this leads to the interaction probability 



(6-15) 



which for the usual AdSs case leads to a scaling with gg and X as a ~ gy^N~^/^. This is a 
rather peculiar prediction for the interaction strength of the string endpoint with the quasi- 
particles. Since the quasi-particles are made out of strings (or branes) some gs dependence is 
to be expected, and the interaction strength indeed vanishes as gg ^ 0. 

In deriving ( |6.15| ) we have assumed that the mean free path time on the stretched horizon 
is identical to that on the boundary. This however, is unlikely to pertain as we explain now. In 
fact, we will argue that on general grounds we should expect that t^fp ~ 1/T on any stretched 
horizon. The logic relies on using dimensional analysis coupled with thermal physics of black 
holes. Generically we expect, 

t.fp = ^^(T/.,^,0 (6.16) 

where ^ is a function of the dimensionless ratios of the length scales available. We have 
fixed the overall normalization to be determined by the thermal scale on physical grounds. 
Furthermore, using the facts that: (i) the dynamics of string probes generically are unaware 
of the Planck scale (to determine which we could for instance use D-brane probes |5^) and, 
(ii) the geometry near a black hole horizon is the Rindler spacetime, which is insensitive to 
the cosmological constant, we can argue that ^ ~ 1, i.e., it is independent of the hierarchy 
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between the various length scales in the problem. More precisely, in the near horizon region, 
<^ 1, the AdSd black hole metric ( p.24 ) reduces to the Rindler metric: 



dsl ^ -f^dP + rff^ + rfXj_2, (6.17) 

where 



i=27rTt, f=^/'—pi, Xd_^ = !j.x^_^_ (6.18) 



Because the metric ( |6.17| ) does not contain any scale such as i, the dynamics of a fundamental 



string in the near horizon region can only depend on Ig (dependence on lp(^d) is excluded as 
in (i)). Therefore, the mean free path time tmfp determined from the dynamics of the string 
can only depend on Ig. However, because tmfp is dimensionless, it should be that tmfp ~ 1, 
which means tmfp ~ 1/T. One can give a more concrete argument by using the argument in 
Appendix applied to the near- horizon geometry ( |6.17|) . 



Now using tjnfp ~ 1/T we can conclude that the interaction probability only depends on 
the ratio of the rf-dimensional Planck scale and string scale: 

which suggests a universal dynamics of the stretched horizon independent of the asymptotics 
of the spacetime. Nevertheless ( |6.19| ) leads to an interaction probability which is a non-trivial 
function of gs as depends non-trivially on the details of the compactification.0 

Clearly, it would be interesting to explore this line of thought further and, for example, 
also find an interpretation for the collision time. However, many of the assumptions we made 
are highly questionable. For example, we ignored backreaction, and only used the quadratic 
part of the Nambu-Goto action. The latter approximation certainly breaks down once we are 
a proper distance ~ Ig away from the horizon. It is also unclear to what extent we can really 
think of the stretched horizon as a gas of almost free quasi-particles. We leave an exploration 
of these issues to future work. 

7 Discussion 

In this paper, we discussed Brownian motion in the holographic context, in order to shed 
light on near-equilibrium dynamics of strongly coupled thermal gauge theories. A useful 
probe exhibiting Brownian motion consists of a fundamental string stretching between the 



■^^ Actually, there is no consensus on where to place the stretched horizon. For example, Refs. [p3|,|54| 
explained some thermodynamical properties of black holes by postulating the existence of quasi-particles 
living on a stretched horizon a distance lp(d) away from the horizon, instead of Ig. In this case with a stretched 
horizon located lp(^d-j away from the horizon, we obtain a simpler result cr ~ 1 instead of ( |6.19 ). This simple 



form of a is appealing, but we do not know of a physical reason to choose one over the other. 
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boundary and the horizon and being randomly excited by the black hole environment. We 
established the relation between the observables associated with such Brownian particle in 
the boundary theory and those of the transverse mode excitations of the fundamental string. 
At the semiclassical level, the modes on the string are thermally excited due to Hawking 
radiation and, consequently, the motion of the boundary Brownian particle is described by a 
Langevin equation, which involves stochastic force and friction. In the bulk, the stochastic 
force corresponds to the random excitation of the string by the Hawking radiation, while 
the friction corresponds to the fact that the excitations on the string get dissipated into the 
horizon. 

Although in this paper we focused on the relation at the semiclassical level between the 
boundary Brownian motion and the dynamics of the fundamental string in the bulk, the 



boundary-bulk dictionary we wrote down in subsection in principle allows one to predict 
the precise correlations of the Hawking radiation quanta beyond the semiclassical approx- 
imation, in terms of the precise correlation functions for the boundary Brownian particle. 
Obtaining the latter of course requires one to compute correlation functions in strongly cou- 
pled plasmas, which is a difficult task. Nevertheless, such a dictionary is an important step 
toward understanding the microphysics underlying the fluid-gravity correspondence. 

One of the particularly interesting results of the current paper is the estimate in subsection 
|3.3| for the time scales associated with the Brownian particle immersed in a CFT plasma: 

Note that setting m = T in troiax gives tmfp, which is a consistency check because a fluid 
particle can be thought of as a Brownian particle with mass ~ T. The fact that tcoii ^ ^^mfp at 
strong coupling A ^ 1 implies that a Brownian particle interacts with many plasma particles 
simultaneously. Because of this, a Brownian particle with mass m <^ \/~\T can thermalize 
in a time much shorter than tcoih the time elapsed in a single process of collision. This is 



reminiscent of the recent conjecture |5^,|5^ that black holes can scramble information very 
fast, whose dual picture is that a degree of freedom in the boundary theory interacts with 
a huge number of other degrees of freedom simultaneously. It would be interesting to study 
this possible connection further. 

Historically, the main achievement of the theory of Brownian motion was the determination 
of the value of the Avogadro constant Ma = 6 x 10^'^ mol~^, which is huge but finite. If Ma were 
infinite, the diffusion constant would be zero and we would not be able to observe Brownian 
motion. The fact that we can observe it in nature gives evidence that Ma is finite and fiuids 
are not continuous but made of molecules. Then, what is the analogue of the Avogadro 
constant in the Brownian motion in the AdS/CFT context we studied, and what is the bulk 
significance of it? In the case of AdS5/SYM4, the macroscopic energy density of the plasma 
scales as = 0{N'^), while the energy carried by a microscopic quantum is of the order of 
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the temperature T = 0{N^). What corresponds to the Avogadro constant is the ratio of 
these, N'^/N^ = N"^. The finiteness of A/a corresponds to the finiteness of A^. In the bulk, on 
the other hand, what corresponds to E is the mass of the black hole, M ~ Rg/G^ ~ 0{G^) 
with Rg the Schwarzschild radius, while T is the Hawking temperature ~ 0{G%). The 
ratio is M /Th = 0{G'j^) = 0{N'^). So, in the bulk, the finiteness of A/a corresponds to the 
finiteness of Gat, or to the fact that the energy carried by a Hawking radiation quantum is 
finite although it is much smaller than the mass of the black hole. 

We considered non-relativistic Brownian motion in the current paper, which is the result 
of the quadratic approximation we made in ( 2.27] ) to the Nambu-Goto action. It would be 
interesting to generalize our treatment to the relativistic case, where Brownian motion and 
its Langevin dynamics are not very well understood; for a recent discussion, see e.g. [^. 
Such a generalization can also be regarded as a generalization of the drag force computations 
of [|0,|22|, which are relativistic because the full Nambu-Goto action was taken into account, to 
non-stationary {u ^ 0) solutions. Also, with such a relativistic formalism, one can presumably 
give a more rigorous derivation of t^fp than the one we did in Appendix |D|. 

As explained above, the stochastic force appearing in the Langevin equation is related to 
the friction term via the fluctuation-dissipation theorem. In the bulk the latter is mimicked 
by the dissipative nature of the event horizon, which is present for all black holes. On the 
other hand, the stochastic term arises due to the Hawking temperature of the black hole; yet 
only non-extremal black holes have finite Bekenstein-Hawking temperature. This leads to the 
naive puzzle that whereas the dissipation is always present, fluctuation is seemingly absent for 
extremal black holes since these have zero temperature. It would be interesting to see whether 
the quantum fluctuations which are present even at zero temperature suffice to account for the 
origin of the stochastic processes. Note that this is not a-priori unreasonable in the AdS / CFT 
context; although quantum processes are 1/A^ suppressed in the large A^ field theory, we had 
to account for the semi-classical Hawking radiation phenomena to see the origin of the random 
force in the Langevin equation. Furthermore, extremal black holes could also be subject to 
super-radiant type instabilities which can effectively mimic the physics of Hawking radiation. 
In fact, this feature has been exploited recently to show how the microstate 'geometries' can 
reproduce some features of the thermal Hawking spectrum ^7}. 

The stochastic random force which drives the long time diffusive motion has a characteris- 
tic dependence on the temperature, which we derived assuming that the system was thermo- 
dynamically stable. As is well known, considering the global as opposed to the Poincare patch 
of AdS provides two distinct black hole solutions at the same temperature — the small black 
hole which has negative specific heat and a large black hole which is in thermal equilibrium 
with the Hawking radiation. To be able to access both these solutions simultaneously one has 
to work in the global AdS geometry which has a compact spatial boundary. The physics of a 
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probe string endpoint in the small black hole background should exhibit marked differences 
from the Brownian motion discussed above, despite the system experiencing the same tem- 
perature. In a finite volume system we naively expect the Brownian process to saturate after 
the time scale t = for in this time the particle has diffused throughout the system. In 

the bulk this presumably corresponds to the string diffusing out completely on the stretched 
horizon and becoming indistinguishable from the thermal atmosphere. This can in fact be 
used to probe the difference between the large and the small black hole. Imagine we normalize 
the physics of the string endpoint on the boundary to correspond to the Brownian motion 
undertaken in the large black hole. Using this as the UV boundary condition for the probe 
string in the small black hole background, we can examine the dynamics of the endpoint at the 
IR stretched horizon. A plausible conjecture for this dynamics is that the fluctuations of the 
string are macroscopically large on the stretched horizon, in fact will have a scale comparable 
to the black hole itself. 
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A Normalized basis 

In this appendix, we discuss the quantization of the action ( p.27|) obtained from the Nambu- 
Goto action and derive the normalized basis of solutions ( p.48| ) to the equation of motion. 

A.l Canonical commutation relations and normalized basis 

The canonical commutation relations for the theory ( |2.27| ): 

SSI = -^J d'^V^)ri^)Gu{x)^^. = t,r, (A.l) 
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are given by 



.2W 



y h \^-^) 

Here, E is a Cauchy surface in the x^ = t^r part of the spacetime ( p.26|) , hij is the metric on S 
induced from (7^1,, and is the future-pointing unit normal to S. For functions f\x),g^{x) 
satisfying the equation of motion (|2.28|) , we can define the following inner product: 

i 



{f,9h 



2na' 



- / dxVhn^Gjj{fd,g'*-d,fg'*). 



(A.3) 



It can be shown that this inner product is independent of the choice of S, just as the standard 
Klein-Gordon inner product This inner product satisfies 



{f,9r = -{r,9l = {9j), (A.4) 
(a/i + 6/2, g)* = a{f,, g) + 6(/2, g), (/, ag^ + hg^Y = a*{f, g,) + b*{f, g^). (A.5) 

It is not difficult to show that the canonical commutation relations ([A.2|) are equivalent to 

(A.6) 



[(/, (5',X)s]s = (/, fi'*)E V/, (7 satisfying the equation of motion (|2.28 



Let {ul^{x)} be a basis of normalized functions satisfying the equation of motion (|2.28|) such 
that 



s. 



a (3, 



0, 



and expand X^ as 



X\x) = ^ [aaul{x) + alul^{xy] . 



Then one can readily show that the condition (|A.6| ) implies 

[tta, ap] = [ai, a^] = 0, [a„, a^] = d^p. 

A. 2 Normalized basis for AdS^ 



(A.7) 



(A.8) 



(A.9) 



As shown in the main text, in the AdSa case, the solution to the equation of motion can be 
written as (see Eq. (p.40|) ) 



-iujt 



(A.IO) 



where B satisfies boundary conditions ( p.41| ) at p = Pc and ( p.44| ) at p = 1 + 2e. The inner 
product ( [A .31 ) for this solution is 



a'p 



2p ^, fp-1 
+ log ^— 

Vp + 1 



1 + p2l/2 



(A.ll) 
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From this, one obtains the normahzed basis: 



2£2culog(l/e) 



fL'-\p) + Bf!:rKp) 



-iujt 



(A.12) 



B Evaluation of displacement squared s^^J^t) 

In this appendix, we evaluate the displacement squared ( ^.5| ), which can be written as: 



where 



h 



dx 



sm 



2 kx 



x{l + a?x'^) — 1 



da; 



1 - e 



ikx 



dx sin^ ^ 



/o X — 1 
and we defined new variables by 

X = 2ttv, 



dxl — e 
\x\ el^l - 1 



P 



\x\{l + a'^x^) el^l - 1 

ikx 



(B.2) 



(B.3) 



27r' /3 

The integrals (p.2| ) can be evaluated using the standard method of deforming the contour 
on the complex x plane. For that, one first replaces |x| with yx^ + e^ with e a small positive 
number, li k > 0, one can then deform the contour to run vertically around the branch cut 
between ie and oo. The resulting integral is simpler than (p.2| ) and, after taking e — > 0, can 
be analytically evaluated. One should also take into account the contribution from the poles 
of the integrand on the imaginary axis. The final result is 



^•4 



eV«Ei(--) + e-*^/"Ei(-) 
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+ ^(1 
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1 . p-27r|fc| 
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2Tia 27ra 

) , 2Fi(i,i-^;2- 1 



2-Ka 



-27r|fc| 
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2™ 



-|fc|/aN 



cot 



2a 



log 



2a sinh vr/c 



(B.4) 



log 



sinh Tik 
irk 



where Ei(2) is the exponential integral, 2-^i(«, /5; 7; -z) is the hypergeometric function, and 
ip{z) = {d/dz)\ogr{z) is the digamma function. For Ei(2;), we take a branch where both 
Ei(x > 0) and Ei(x < 0) are real. 



41 



If Pc ^ 1 and thus a ^ 1, one can use the expressions ( [B.4|) to derive the following 

lavinr- 



behavior: 



7ik + 0{\ogk) 



C(aO) 

Trk + 0{\ogk) 



{k < a) 
{k > a) 



(B.5) 



Therefore, if pc ^ 1, s'^eg^it) has the following behavior: 



{t « /?) 
» P) 



(B.6) 



C Distribution of momentum p 



In this appendix, we compute the probability distribution of the momentum p = mx, where 
X is the position of the string endpoint at the UV cut-off p = pc, and show that it is exactly 
equal to the Maxwell-Boltzmann distribution.^ 
From (|2.51|) , the momentum of the particle is 



p = mx{t) 



im 



E 



log(l/e) 



1 — iu ( p, 



1 - ipcV \pc + 1 



iu/2 



-iujt 



a„, — h.c. 



(c.i: 



We would like to know the probability distribution f[p) of p. By definition, 

/oo 
dpe'^^f{p). 
■oo 

Namely, f{p) is the Fourier transform of (e*^^). So, what we want to compute is 

l-iu f pc- 



(C.2) 



(:e-^:) = (:exp<(^5: 

w>0 



2a'(3uo 
log(l/e) 



1 - ipcV \pc + 1 



iv/2 



-iuit 



— h.c. 



(C.3) 



where we regularized the operator by normal ordering. The expectation value is with respect 
to the density matrix ( p.52|) . Using the identity 



tr 



1 - e-P^ 



exp(- 



we can compute 



(:e^P«:) =Cexp 



2^/^2 /-oo 



duu{l + 



{l + plu^W-^ 



(C.4) 



(C.5) 



^^Here, we will ignore the fact that the mass of the quark gets corrected in thermal medium and make 
a crude estimate by using the bare mass ( 2.39| ). 
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where C is a constant independent of ^ and we rewrote the sum over uj by an integral using 
( p.45 ). This integral can be evaluated by deforming the contour in the complex plane, just 
as we did for ( |B.2|) , the result being 



dvviX + 



(l + p2i.2)(e2.._i) 



TTCOt ^ ) - 

Pc Pc 



Pc 



21ogpc 



+ 



24p2 



(C.6) 



where ip{z) = (d/dz) log r{z) is the digamma function. Using this expression, it is easy to 
show that, for large pc, 

duu{i + u^) _ 1 

X (l+p,V2)(e2--l) ~4ic^ 
Therefore, from ( |C.5| ), we obtain 

for large pc, where we used ( p.39| ). By Fourier transforming, 
fip)- 



(C.7) 



(C.8) 



2ii 



Er, 



2m 



(C.9) 



This is exactly the Maxwell-Boltzmann distribution of particles with energy Ep. Therefore, 
for large pc, the endpoint of the string behaves like a non-relativistic particle with mass m 
immersed in a thermal bath of temperature T. 



D Mean free path time t^ifp 



In subsection |3.3| , we discussed the time scales associated with Brownian motion: the relax- 
ation time trciax and the collision duration time tcou- In this Appendix, we evaluate tmfp, the 
mean free path time, or the typical time between two collisions, using the correlators of the 
random force R{t) in the case of AdSs. We argue that the mean free path time is given by 

(D.l) 

where A ~ /a'"^ is the 't Hooft coupling, although we are unable to give a rigorous derivation. 
We expect that this holds in more general cases, including AdSs. 



In subsection p.l| , we discuss how to determine characteristic time scales from correlators 
in general. In subsection p.2| , we compute tmfp for the Brownian motion in the case of AdSa. 



43 



D.l Correlators and time scales 



Consider a stochastic quantity R{t) whose functional form consists of many pulses randomly 
distributed. Let the form of a single pulse be f(t), with width A and amplitude A. Fur- 
thermore, assume that the pulses come with random signs. If we have k pulses at t = 
{i = 1,2, . . . , k), then R{t) is given by 

k 

R{t) = Y,^,f{t-U), (D.2) 

i=l 

where ej = ±1 are random signs. For a schematic picture, see Figure |^. Let us assume that 

R(t) 

Figure 5: A sample of the stochastic variable R{t), which consists of many pulses 
randomly distributed. 



the distribution of pulses obeys the Poisson distribution. Namely, the probability that there 
are k pulses in an interval of length r, say [0, r], is given by 

Pk{r) = e-'^^-^. (D.3) 

Here, /i is the number of pulses per unit time. In other words, is the average distance 
between two pulses. We do not assume that the pulses are well separated; namely, we do 
not assume A ^ Later, we will identify R{t) with the random force in the Langevin 

equation; the pulses are contributions from a collision with a fluid particle, and therefore 

The 2-point function for R can be written as 

oo ( \k ^ 

{R{t)Rit')) = E ^ kT ^ ^''''^^^ ~ " 

fc=l i,j=^ 

where we assumed t,t' G [0, r] and ( )^ is the statistical average when there are k pulses 
during [0,r]. Because k pulses are randomly and independently distributed in the interval 
[0, r], this expectation value is computed as 



k 

" k k 

dti ■ ■ ■ dt, 



. (D.5) 
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Here, the second term vanishes because {^i^j)f^ = for i 7^ j. Therefore, one readily computes 
{^iejfit - U)f{t' - t,)), = - r dtJit - U)f{t' - h) 

A „-i "''0 

/OO 7. 
dh fit -t'- h)f{-h) = -F{t - t'). (D.6) 



k 



Here, in going to the second hne, we took r to be much larger than the width A of /(t), 
which is always possible because r is arbitrary. Substituting this back into (|D.4| ), we find 



{R{t)R{t')) = fiF{t-t'). (D.7) 
In a similar way, one can compute the following 4-point function: 

00 f \k ^ 

k=l i,j,m,n=l 

(D.8) 

Again, the expectation value {eiejem^n))^ vanishes unless some of m, n are equal. The 
possibilities are i = j^m = n, i = my^j = n, i = nj^j = m, and i = j = m = n. 
Therefore, 

k 

Y {^i^j^menf{t - U)f{t - t^)f{t' - tjf{t' - tn))^ 
i,j,m,n=l 
k 

[Y [fit - t'd'fit' - + 2f{t - u)f{t' - u)f{t - t,)f{t' - 1,)] 
k 



r2 



i=l 

r>00 



00 



dhdh [f{t - hffit' - hf + fit - h)f{t' - h)f{t - t2)f{t' - t2)] 

k f°° 

+ -/ dt^f{t-tiff{t'-t,)\ (D.9) 

T J -00 

Substituting this back into ( p.8|) , we obtain 



{R\t)R^{t')) = /i'[F(0)' + Fit - t'f] +12 dufit-t' - uff{-uf 

J ~oo 

/OO 
duf{t-t'-uff{-uf. (D.IO) 
-00 

For example, consider the following shape function for the pulse 
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Then, one computes 

{R{t)R{0)) = v^A/iAV*'/^^', (D.12) 

{R'm'm = {R\t)){R'{0)) + {R{t)R{0)f + -=l—{R{t)R{0)f. (D.13) 

yiTT A /i 

Therefore, if we know the behavior of {R{t)R{0)), we can read off A and fiA from (p.l2|) . 
If we further know (i?^(t)i?^(0)) then, from (|D.13| ), we can read off fi. In particular, if we 
denote the last term in ( p.l3|) by {R'^)', then 

^-^^!i^A. (D.14) 



This result ( p.l4D is expected to be true for other forms of f{t), not just for the Gaussian 
case ( pup . 

Note that the treatment above is classical. If R{t) is a quantum operator, we should con- 
sider the classical part of the correlators by appropriately subtracting quantum divergences. 

D.2 Evaluation of tnifp for Brownian motion 

Let us evaluate t^fp for the boundary Brownian motion, by identifying the stochastic function 
R{t) in the previous subsection with the random force appearing in the Langevin equation. 
A pulse f{t) corresponds to the contribution from a collision with a single plasma particle. A 
is the time elapsed in a single collision, namely A = tcou, while l//i is the time between two 
collisions, namely = tmfp- 

Using Eq. ( p.51| ) as well as the relations p = mx and R{uj) = p{uj)/fi{uj), where fi{uj) is 
given by ( |3.23|) , we can write the random force R{t) as 



^ \l a'(3-^\og{l/e) 1 - tu/pc \pc + l ' 



uj>0 



Because free harmonic oscillators, it is easy to show that 

{■.R\t)R\0):) = (:i?2(t):)(:i?2(o):) + (; i?(t)i?(0) :)l (D.16) 

Here, we are considering the normal-ordered correlators because the result of the previous 
subsection applies to the classical piece of correlators; henceforth, normal ordering of operators 
will be understood. By comparing (P.16| ) with ( p.l3|) , we appear to have tmfp = 0. However, 



this is due to the non-relativistic approximation we made in ( |2.27|) when we expanded the 
Nambu-Goto action up to quadratic order. If we keep the next order (quartic) terms, we 
obtain the following additional contribution to the Hamiltonian: 
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H^^^ = — I dr 



\dtXf r^h{r 
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h{r) £4 



idrX) 



Mr) = !-(-) , (D.17) 
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where = (1 + 2e)rH, e <^ 1. This corresponds to the first relativistic correction to the 
non-relativistic action S^q. In the presence of this interaction, there is an extra contribution 
to the correlator {R'^{t)R'^{Q)) coming; from the contractions with the terms in H^^\ If we 
consider the case with t = 0, we have 



0' 



(D.18) 



where ( )q is the expectation value with respect to the quadratic action 5'^;^^, i.e., it is the 
expectation value with respect to the density matrix (|2.52|) . 

So, let us evaluate the last term in (|D.18| ), which will be denoted by {R'^)' ■ Using the 
expansions ( |2.47|) and (|D.15| ), the explicit expression for (i?^)' is 



16na' 



\iiJl,...,UJ4>0 



r'^h 



£4 



. (D.19) 



~oj>0 ~ ' uj>0 ' J / 

There are many terms coming from the expansion of this. Let us focus on the following term 
in particular: 



UJ1,...,Ll!4 UJ ,...,Ll1' 



■'Is 



"^1 "-^2 "''^3 

(D.20) 

There are various ways to contract Let us take the term obtained by contracting a^^ 

against a}^,, or against a^^', where z = 1, . . . ,4. Other contractions give similar contribu- 
tions. This particular contraction gives the following: 
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(D.21) 



where the Bose-Einstein factor l/(e^'^' — 1) has effectively cut off the Ui sum at f3~^. From 
now on, we do not keep track of numerical factors. 

We would like to evaluate the r integral in (p.21|) . Because the integrand in ( p.21|) has a 
second order pole aX r = th due to the dominant contribution comes from r ~ ~ th- 
For a while, let us instead consider the case where there is a first order pole at r = rn, by 
replacing h~'^ by . From ( |2.50| ) and (|2.31|) , near r = r^^. 



a'(3 



2£2a;log(l/e) 



dr r% , 
— ^ —^nr 
h & * 
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2n2 



-dr^ 



(D.22) 
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Therefore, the integral in ( p.21|) (with h ^ replaced hj h ^) is 



a 



/2 



£2[log(l/e)]V^1^2W3W4 7_Alog(i) 

^ j-giojar, _|_ gi9ij3g-iaJ3r*j |-g-iaJ4r. _j_ ^-461^4 gia;4r*j +0)2+1^3-^^4)* 

The dominant part in the e — limit can be easily evaluated by noting that 
/ dr, e'"^* = 5^ ^ log f -) + (finite as e ^ 0) . 

For example, by collecting the first terms in the four pairs of the brackets in ( p.23| ), one finds 



(D.23) 



(D.24) 



£2 log(l/e) y/UiUJ2i03UJ4 



(D.25) 



Note that the finite part in ( p.24| ) does not survive in the e — limit. If we plug this result 
back into (p^) , using the explicit expression for in (|D.15| ), we find 



i</3- 



m/Pc \Pc + l 



Pc 



iui/2- 



(D.26) 



Here, "[2]" denotes the previous factor with z/i replaced by i^2- "[3]" and "[1 + 2 + 3]" are 
similar. By rewriting the sum by integral using ( p.45| ) and using the fact that 1, this is 
estimated as 

/C2 /■ ^2 1 



duoi duo2 duJ3 



(D.27) 



where we used the relation A ~ ^^/a'^. There are many other terms we did not discuss, such 
as other contractions of (p.20|) , but these will not affect this estimate. 

However, of course, this is not precisely what we wanted to evaluate; we have replaced 
in ( p.21|) by h~^. However, using instead will change the above discussion, because the 
r integral around r = th will now give a power (~ 1/e) divergence instead of the logarithmic 
divergence we had in ( p.24|) . This log divergence was important in obtaining the result 
( p.27|) , because this log divergence was canceled against the normalization factor in ~ 
[log(l/e)]~^/^. What we have forgotten is that, if we include the quartic correction H^^\ we 
should also consider corrections to the normalized basis u^, which presumably introduces 
a normalization factor that goes as e^^^. This corrected normalization factor should cancel 
against the power divergence coming from thus giving a finite result, which should give 
( p.27| ) at the end of the day — namely, 

1 



(D.28) 
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Whatever the modifications due to the quartic term are, the dominant contribution comes 
from the region r ~ r^/ and quantities such as Tc or m cannot enter the final result. Also, the 
relativistic correction must come with a factor of a' ~ There being no other available 

quantities, (-R^)' must be proportional to (|D.2^ ). A fully relativistic formalism in which one 
can rigorously and explicitly show ( p.28|) is beyond the scope of the current paper. We leave 
development of such a formalism for future research. 
From ( 3.35|) , we have 

(fi^>-°(*-0)^^^^. (D.29) 



Therefore, using the formula ( p.l4|) with A = /3, we obtain 



It is satisfactory that this does not depend on the properties of the Brownian particle probe 
such as m, because tmfp is a time scale associated with the fluid itself. 

E Solving equation of motion for general d using match- 
ing technique 



In this appendix, we solve the wave equation ( |5.6| ) for general dimensions using the matching 
technique for low frequencies u <^ T. We would like to obtain a solution which is purely 
ingoing at the horizon and, in particular, determine its behavior near the boundary. 

Let us denote by X~{r) this particular solution of the wave equation (|5.6|) which obeys 
the purely ingoing boundary condition at the horizon r = th- To determine it, let us consider 
three regions: (A) a near horizon region with r Th and V{r) ^ cu^, (B) an intermediate 
region with V{r) ^ uj"^, and (C) an asymptotic region with r ^ rn- The idea is to consider 
the approximate solutions in each of the three regions, and to match these to each other. For 



more details, see 15811 and references therein. As before, we define 



P = — y = — • E.i 

In terms of these parameters the constraints on the different regions under consideration, 
V{r) <^ uj"^ and V{r) ^ u;^ respectively translate to p — 1 ^ and p — 1 ^ z/^. Furthermore, 
the low frequency condition, u; ^ T, can be written as z/ ^ 1. 

In region A, where p — 1 <^ z^^ and we can drop the potential V{r) from ( p.6| ), the linearly 
independent solutions are 



X^{r) =e±^^"* ~exp 



:^log(p-l) 



(E.2) 
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The purely ingoing solution is X^{r) = e"*"^^*. Now, since ^ 1 we can focus on a region 
slightly away from the horizon (still remaining in region A), such that exp(— ^ p — 1 ^ 
z/^. Here we can approximate the purely ingoing solution as 

i V 



d 



log(p - 1). 



(E.3) 



In the asymptotic region C, where p ^ 1, we can approximate h{r) ~ 1. The linearly 
independent solutions of (^) are then 



The general solution can be written as 

Xc 

which, for p 1, can be expanded as 

Xc = (C+ 



It 



p 



+ C X(j , 



c-)|i + ^ + ... 



3p3 



(E.4) 



(E.5) 



(E.6) 



Finally, in region B, where p—l^v^ and we can drop cu^ from ( |5.6| ), leading then to the 
general solution 

dp' 



xb{p)=b,+b, r 

J oo 



where Bi and B2 are two integration constants. For p ~ 1 (but still p — 1 ^ i/^), we can 
approximate h{p) ~ (rf — l)(p — 1) and ( |E.7| ) gives 

1 



d 



log(p 



(E.8) 



where 6 is a constant independent of v whose precise value is not relevant for our purpose. 

We now have the solutions in the three regions A-C; by matching them across the domains 
of overlap we can relate the various constants of integration. To begin with we determine Bi 
and B2 by matching ( [E.8|) in region B with the solution ( [E.3|) in region A, obtaining 



Bi 



i h u. 



Bo = -iu. 



(E.9) 



To determine we expend the solution in region B ( |E.7| ) for p ^ 1 and match it to that in 
region C ( |E.6| ) leading to 



B^=C+ + C- 



Bo 



(E.IO) 



It must be borne in mind that we have performed the matching only in the small frequency 
limit -C 1 and as a result should trust the expressions only at the leading order in u. Solving 
( [E.9|) and ( PH.IOI) , we finally find that the purely ingoing solution behaves at large p as 



Xj(p)=C+X+ + C"Xj, 



- {l± — + ibu 



(E.ll) 
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